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> Abstract. We present a new prescription for analysing cosmological perturbations in a more- 

general class of scalar-field dark-energy models where the energy-momentum tensor has an 
imperfect-fluid form. This class includes Brans-Dicke models, f(R) gravity, theories with 
t-h kinetic gravity braiding and generalised galileons. We employ the intuitive language of flu- 

ids, allowing us to explicitly maintain a dependence on physical and potentially measurable 
properties. We demonstrate that hydrodynamics is not always a valid description for de- 
scribing cosmological perturbations in general scalar-field theories and present a consistent 
alternative that nonetheless utilises the fluid language. 

We apply this approach explicitly to a worked example: k-essence non-minimally coupled 
to gravity. This is the simplest case which captures the essential new features of these 
imperfect-fluid models. We demonstrate the generic existence of a new scale separating re- 
gimes where the fluid is perfect and imperfect. We obtain the equations for the evolution 
of dark-energy density perturbations in both these regimes. The model also features two 
other known scales: the Compton scale related to the breaking of shift symmetry and the 
Jeans scale which we show is determined by the speed of propagation of small scalar-field 
perturbations, i.e. causality, as opposed to the frequently used definition of the ratio of the 
pressure and energy-density perturbations. 
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1 Introduction 

Many models, better or worse motivated by fundamental physics, have been proposed as 
an alternative to the cosmological constant for providing a mechanism to accelerate the 
expansion of the universe at late times. Perhaps the simplest, and therefore most studied, 
are models which feature a single extra scalar degree of freedom, which is treated as classical. 
The archetypal example of this dark energy ("DE") are quintessence models [1, 2], where 
a canonical scalar field rolls in a potential resulting a time-varying equation of state and 
allowing for (small) inhomogeneities in the dark-energy fluid. On the other hand, f{R) 
gravity models [3], where the Einstein-Hilbert action for gravity is modified, allow for a 
different response of the metric to the presence of matter from that in general relativity 
("GR") and usually are discussed as examples of modified gravity. However, these models 
can also be reformulated as a particular subclass of Brans-Dicke theories featuring a single 
extra scalar, albeit non-minimally coupled to gravity [4, 5]. It was quickly understood that 
whereas the background expansion history can be indistinguishable between various classes of 
these models, the evolution of dark-matter perturbations can differ significantly and therefore 
formation of large-scale structure provides the key to understanding the role played by these 
theories, if any, in accelerating the universe [6]. 

k- essence was introduced as an extension of quintessence models to non-canonical kinetic 
terms [7-9] and can be interpreted as a perfect fluid with a speed of sound different from that 
of light [10] . This allows for a Jeans scale smaller than the cosmological horizon and therefore 
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for DE which can cluster significantly. In particular, the discovery of generic attractors where 
the equation of state approaches that of the vacuum while the sound speed tends to zero [11] 
also spurred the discussion of such models as describing both dark energy and dark matter 
through a single degree of freedom [12-14]. One tends to think of these dark-energy models 
as perfect fluids, obeying hydrodynamics and studies the evolution of linear perturbations 
by considering the conservation of the perturbed energy-momentum tensor [15-17]. 

On the other hand, when studying large-scale structure in f{R) theories [18, 19], one 
tends to think in terms of a modification of the response of the two scalar gravitational 
potentials, <]? and fy, to matter perturbations. In particular the generic result for this class 
of theories is that the lensing potential, <£ — *!', is "unaffected" and depends only on the 
matter perturbations and an evolving Newton's constant, while the modification of gravity 
creates anisotropic stress giving 2$ + \& = 1 . These relationships are simple and therefore 
are frequently generalised in order to constrain such theories from data. A large number of 
phenomenological parameterisations of models of modified gravity exist [6, 21-28] and usually 
focus on directly modifying the relationships between the potentials and matter perturbations 
without necessarily asking what this implies for the underlying model. 2 Considering a lack 
of well-motivated models, one may argue that this is not a significant problem. 

f(R) models and this sort of modifications of gravity are quite constrained by Solar- 
System tests as a result of their being equivalent to fifth forces. One can find viable models 
[30-32] which can evade the constraints through the chameleon effect [33, 34], or one can pro- 
pose that the extra scalar does not couple to baryons, but only to dark matter [35, 36]. These 
types of interacting DE models can have dynamics on large scales very similar to modified 
gravity models, since the contribution of baryons to the energy density is very subdominant. 
However, philosophically they are no longer "modified gravity" since the coupling is not uni- 
versal — it no longer couples purely to the energy-momentum tensor. On the other hand, the 
full breadth of such models is much wider than for f{R) theories, which are but a subclass. 
When studying interacting dark energy, one tends to treat the scalar field differently from 
matter: the equation of motion for the scalar is usually parameterised and solved, while the 
dark matter and baryons are treated as a hydrodynamical fluid. 

Another class of scalar-field models was discovered by considering gravity modifica- 
tions arising from a potential higher-dimensional embedding of our universe [37-39]. In the 
decoupling limit [40-42], this DGP gravity reduces to a scalar theory, where the scalar Lag- 
rangian contains second derivatives of the scalar field and therefore is not in the k- essence 
class. A further generalisation of such scalar-field theories featuring second derivatives was 
then obtained in refs [43, 44] and named the galileon. Such theories are not consistent on 
general curved backgrounds, and the appropriate covariantisation was obtained in [45, 46]. 
Kinetic gravity braiding, a generalisation of the DGP decoupling limit, was independently 
found in refs [47, 48] and soon it was proven that the most general theory of a single scalar 
with no more than second derivatives is a similar generalisation of the galileons [49] . Indeed 
these generalised galileons were discovered a long time ago [50] and largely ignored. The 
higher-dimensional origin of the galileon theories was studied in refs [51-54], while the ga- 
lileon theories also have appeared in the decoupling limit of consistent massive-gravity as 
actions for the helicity-0 mode of the massive graviton [55-57] . 

For an investigation of the relation between anisotropic stress and cosmological evolution in non-linear 
gravity models see ref [20]. 

2 Although see ref. [29] where a very general parameterisation is given based on symmetries of the Einstein 
equations. 
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These theories all contain second derivatives in the action and therefore have energy- 
momentum tensors ("EMT") which do not have the form of a perfect fluid [58, 59] and 
therefore provide for new phenomenology. For example, the null-energy condition can be 
violated stably, resulting in explicitly ghost-free phantom DE models or even permitting the 
construction of initial phases of evolution of the universe alternative to inflation [60-62]. 
Equations for linear cosmological perturbations for the most general of such theories were 
derived in [63]. 

What strikes in the above discussion is the multitude of approaches that are used to 
study perturbations, each chosen depending on the underlying Lagrangian for the scalar. 
However, all the models discussed here are structurally extremely similar. A common frame- 
work for describing perturbations would allow for a single consistent method of constraining 
the class of models of DE comprising a single scalar. One could argue that solving the equa- 
tion of motion for the scalar provides exactly such a framework. However, if we assume that 
the DE cannot be studied directly, but only through its effect on the gravitational field in 
which baryons and light propagate, the equation of motion contains too much information. 
Gravity couples only to the EMT and therefore will mostly depend on the fluid properties of 
the DE and not on the precise value of the scalar field. 

The case of k- essence is a pertinent example. Any particular choice within this class 
of Lagrangians gives a particular expansion history for the cosmological background, with 
a particular evolution of the equation-of-state parameter, w. This solution (based on the 
choice of Lagrangian) then also determines the speed of propagation of the scalar-field per- 
turbations, Cg . Given the history of these two variables, the evolution of linear energy-density 
perturbations of a k-essence DE is completely determined at all scales. Any k-essence model 
which happens to have the same history of evolution of w and will have exactly the same 
evolution of linear perturbations, even if the two Lagrangians are very different and only 
happen to have coincided in this way on two individual solutions. There is therefore no way 
to reconstruct a Lagrangian even limited to this single class of theories, even if we ignore the 
fact that the data that can be obtained will only ever be limited to some maximal precision 
and we will only be able to measure the properties of the DE indirectly. Unless DE interacts 
directly with baryonic matter in a way that is not equivalent to redefining a gravitational 
metric, we will only be able to measure the properties of the DE perturbations indirectly by 
looking at scale-dependent modifications of the DM power spectrum. Given this, the best we 
can hope for is to measure w through its effect on geometry and the sound speed of DE by 
detecting some sort of scale of transition in behaviour that would be consistent with a Jeans 
length [64, 65]. 

A similar situation occurs when dark energy is described by the class of f{R) theories. 
Here, there appears to be only one scale, that of the (time-varying) Compton wavelength of 
the scalar, the chameleon effect allowing for the restoration of GR in dense regions through 
non-linearities notwithstanding. On both sides of this scale, the evolution of linear perturb- 
ations is completely determined by the fact that we have chosen to consider this limited 
class of Lagrangians, just as it was for k-essence. Indeed, specifying the evolution of w and 
of the Compton scale gives us the totality of the information that can be gleaned from the 
background and the linear perturbations. 

In fact, the standard discussion of f{R) gravity neglects two more scales which are 
implicitly contained in the theory. Firstly, there is still the Jeans scale. However, in this 
case, the speed of sound is equal to that of light and therefore the Jeans scale lies at the 
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cosmological horizon and is thus unobservable. Secondly, there exists a new scale controlling 
whether the fluid is perfect of imperfect, which will be the subject of much of this paper. 
f(R) models are a limit where this scale lies at infinity and therefore it is also unobservable 
there. 

The main aim of this paper is to provide a prescription for describing linear perturbations 
which as far as possible depends only on physical properties of the model, such as the equation 
of state or the sound speed, which one could hope to measure. Our approach is to describe the 
DE through the conservation equations for its EMT, rather than the equation of motion. This 
provides the evolution equations for the energy-density perturbations which are then directly 
connected to the gravitational potentials. For any theory, these equations can be solved 
provided we supply two closure relations: relationships of the pressure perturbation and 
anisotropic stress to the energy-density perturbation. The solutions for linear perturbations 
of any fluid are always fully determined by these two functions. 3 The question we answer 
is how to properly obtain these relationships for a model featuring a single scalar degree of 
freedom. 

We show explicitly that even for general k-essence models these closure relations, which 
are well known [66], are not hydro dynamical and one should never a priori assume hydro- 
dynamical behaviour in more general cases. However, on any time slice all components of 
the EMT are determined by the values of the degrees of freedom at that moment in time. 
Thus, for any class of models based on degrees of freedom with classical equations of mo- 
tion, we can obtain such closure relations. Indeed these relations define the totality of the 
properties of the dark energy at the linear level of perturbations. We provide a form of these 
closure relations for any scalar-field theory and we explicitly calculate them for k- essence 
coupled non-minimally to gravity. This particular class contains both k-essence and f(R) 
theories and as a result of the non-minimal coupling to gravity features second-derivatives of 
the scalar in the EMT in the Jordan frame. In this sense, it mimics the form of the EMT 
possessed by more-general galileon theories and therefore is a useful example. 

We show that in general scalar-field models we have three scales: the Jeans scale, the 
Compton scale and a new scale determining whether the fluid is in a perfect or imperfect re- 
gime. This is in addition to the cosmological-horizon scale. These three scales are essentially 
independent and can lie in any order and are determined by the Lagrangian class and the 
particular background solution. This new transition scale appears in all scalar-field models 
more general than k- essence and is determined by the relative importance of the second- 
derivative terms and the k-essence terms in the perturbed EMT. For example, it is only 
inside this scale that the DE fluid can carry anisotropic stress at all. This scale is not visible 
in the equation of motion for the scalar, but only appears when the EMT is considered. 

As an example of the power of our prescription, we obtain the equations for the evolu- 
tion of energy-density perturbations for the non-minimally coupled k-essence DE model and 
solve them analytically under the assumption that the background evolution exhibits scaling 
behaviour (i.e. has a constant equation-of-state parameter). 

The paper is structured as follows: we start in section 2 by describing a formalism 
for decomposing the EMT into the fluid variables covariantly. This description allows us 
to obtain exact forms for, among others, the energy density and pressure in terms of the 

3 Some EMT non-conservation issues related to the evolution of the effective Planck mass notwithstanding. 



-4- 



scalar-field, which can then be perturbed directly. In this way, we obtain the conservation 
equations for the linear perturbations of the EMT in section 2.2. In section 2.3 we discuss 
what hydrodynamics implies for the relationship between pressure and energy density in 
perfect fluids, and by deriving this relationship for general k-essence models in section 2.4 
we show that the hydro dynamical relations are not obeyed. We then propose a form for 
these closure relations which would be obeyed by any DE model comprising a single scalar 
field in section 2.5. In section 3, we turn to a non-minimally coupled k-essence as worked 
example. We discuss the general properties of this theory, before turning to calculate the 
closure relations between the pressure, energy density and anisotropy perturbations in section 
3.3. We demonstrate explicitly how to calculate them and show that they are of the form 
contained within the general parameterisation introduced earlier. We explicitly demonstrate 
the existence of the three scales in the problem and derive the equation for the evolution of 
the density contrast in these models. We summarise our findings in section 4. 



2 Fluid formalism 



It is well known that theories such as k- essence can be much more intuitively described in 
terms of relativistic hydrodynamics [7]. We are going to employ this language extensively 
and show that it provides a very natural framework in which to discuss theories with con- 
tributions to the energy-momentum tensor containing second derivatives. This language was 
also previously employed to study scalar field theories with kinetic gravity braiding [58] (to 
where we refer the reader for a more detailed exposition of this formalism), but was developed 
for relativistic potential flows much earlier [67-70]. We will proceed in this manner, leaving 
aside for the moment the question of whether the more general scalar-field theories we would 
like to describe are in fact hydrodynamical fluids. 

2.1 Standard definitions 

In this discussion, we are going to be describing a scalar-field theory. Following standard 
practice, we will define a canonical kinetic term for the scalar, 4 

X = ~d^d»4>. (2.1) 

In cosmology one usually takes gradients as being time-like. It is only with this choice that 
the fluid description can be valid, thus we will always assume that X > 0. The scalar- 
field gradient now gives a time-like vector field, which when appropriately normalised can be 
identified with a velocity field, 

This velocity field is usually referred to as the fluid rest-frame, but in the case of the more 
general models we discuss, this is not always the case: there remains a non- vanishing energy 
flow in the frame moving with the velocity u^. We will therefore call the frame defined by 
(2.2) the scalar frame. We comment further on this frame choice in Appendix B. 

Given the velocity field, we can define a derivative along vP" (a material derivative) 



d 

dr 

4 We employ the metric signature convention ( — h ++ 



u^V M = — , (2.3) 
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thus making r the proper time of an observer comoving with the scalar frame. We will 
assume that this time derivative be positive, and therefore be future directed. This allows 
us to think of eft as a clock. 

We can then proceed and define the transverse projector, 

-L/ii/= 9/iu + Ufj,u u , (2.4) 

which plays the role of the first fundamental form in the hypersurfaces : <p{x) = const, i.e. 
the spatial metric for the scalar-frame observer. The projector (2.4) allows us to decompose 
vectors and gradients into time and spatial parts as observed in the scalar frame, 

V M = -u M « A V A + ±J V A = -u„— + V M . (2.5) 

We will refer to as the spatial derivative. We can then proceed to decompose the gradient 
of the vector field by defining the acceleration, a spatial vector field 



d7 



u M = u A V A ii At , (2.6) 



which can also be expressed in terms of the scalar field 

a M = -i.JV A ln^, (2.7) 
and performing the standard kinematical decomposition 

V M u^ = -u^dv + a^ u + - 9 , (2.8) 
where 9 = V^u^ is the expansion and 

=± x itl V x u v) - -±„ u 9 (2.9) 

is the fully spatial (transverse) and traceless symmetric shear tensor. The rotation tensor 
(the antisymmetric spatial part of (2.8)) vanishes by virtue of the Frobenius theorem, since 
Ufj, is a gradient of a scalar. 

In ref. [58], it was proposed that instead of the canonical scalar kinetic term, X, a more 
physical variable is in fact 

m^^/V^f (2.10) 

On the cosmological background, which will be the focus of our discussion, the variable m 
reduces to the derivative of the scalar with respect to the coordinate time, d^/di. It will 
also be helpful to give the decomposition of the second derivative in terms of the kinematical 
variables: 

dvn ( I \ 

V M V^0 = --^u^Uv - m ia^Uy + u^a u - - - ± fiu 9) . (2.11) 

In the hydrodynamical picture, the variable m can be interpreted as playing the role of 
the chemical potential for the charges which make up the fluid. In this picture, a scalar- 
field theory is a thermodynamical system at zero temperature, with m acting as a chemical 
potential and any dependence on <f> representing an explicit dependence of the properties on 
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the internal clock. As argued in [71-73], scalar-field theories such as k-essence are perfect 
hydro dynamical fluids only when they are shift symmetric, with the action invariant under 
the transformation <j)(x) — > 4>{x) + const, i.e. hydrodynamics is only valid when there is no 
explicit dependence on the clock carried by the fluid. We will return to this point in the 
discussion of the validity of hydrodynamics in sections 2.3 and 2.4. 

The configuration of a minimally coupled k- essence theory is fully specified by the pair 
(<p,m). Theories with second derivatives in the Lagrangian, such as kinetic gravity braiding 
[58] are somewhat more complicated, since their configuration can also depend on quantities 
such as the expansion 9. Typically, the energy momentum tensor of these theories will contain 
second derivatives of the scalar. We will discuss this sort of example in section 3. 

2.2 Evolution equations from EMT conservation 

At this stage, it is useful to jump ahead and be a little more specific. In scalar models coupled 
non-minimally to gravity, the fundamental Planck mass, Mp\ will in general be modified by a 
function of the scalar, let us call it x, through a coupling e H R in the Lagrangian. This would 
mean that the energy momentum tensor ("EMT") for the DE would in general contain a 
term proportional to the Einstein tensor, (1 — e K )G^ u . When the Einstein tensor is solved for 
through the Einstein equations, the energy-momentum tensor for the dark energy depends 
explicitly on the dark-matter configuration, which we would like to think of as a separate 
degree of freedom with its own dynamics (in fact, potentially multiple degrees of freedom: 
dark matter, baryons, radiation, etc.). The prescription we introduce in this paper calls for 
separating out this dependence, so as to leave the DE EMT as much as possible independent 
of the particular configuration of the external matter. We will therefore explicitly solve for 
the Einstein tensor and introduce an EMT, T^ v with the Einstein tensor removed. The 
Einstein equations are then given by 

Gp, = TX+e-*Tfi, (2.12) 

where we have chosen units where Mpj = 87tGn = 1. Tff is the possible EMT external 
to the dark-energy system, for example representing dark matter, which would be explicitly 
conserved, VIp 1 = 0. All the fluid variables relating to the external matter will appear 
suppressed by e~ K '. The price for such a definition is the fact that T^ v is not conserved in 
the presence of matter external to the system, 

V^T* = e-^Vx • (2.13) 

In this paper, will consider the isolated dark-energy system, setting = 0, only making 
comments in passing when the presence of external matter modifies significantly our results. 
We will return to study the combined dark-energy-dark-matter system in another work. 

A natural question to ask is why we are not performing the analysis in the Einstein 
frame, which is the typical way in which non-minimally coupled theories are usually dis- 
cussed [74]. Indeed for the model we present in the worked example in section 3, such a 
transformation can be performed and the analysis of such a setup is standard: in the Ein- 
stein frame we are dealing with a normal perfect fluid. However, more general models of 
DE (kinetic gravity braiding, generalised galileons) feature an intrinsic deviation in the EMT 
from the perfect-fluid form which cannot be removed by a redefinition of frame. The Ein- 
stein frame defined as the frame where the kinetic terms of the metric and of the scalar field 
are unmixed, does not necessarily exist in more general models. Our motivation is that the 
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simple model presented in our worked example, in the Jordan frame, has many features in 
common with the general galileon models and therefore our prescription as employed in the 
Jordan frame will transfer directly to more complex situations. 

In general, the baryons, DM and gravity could all have different frames in which they 
have minimal couplings. We can always pick one of those frames as the basis for the analysis. 
Since most of our observations are to do with baryons, it simplifies matters to assume that 
they move on geodesies of some metric. This is the metric in which they are minimally 
coupled. But in this frame, the DM and gravity could both have some sort potentially dif- 
ferent non-minimal coupling for the scalar. 

Once we have defined a velocity field, we can decompose the EMT into the fluid variables 
that will be seen by the observer O comoving with velocity u 11 : 

T* v = £u^u u + J_ M „ V + Ufj,q v + u u q^ + t^ u , (2.14) 

where £ = T^ u u^u u is the energy density observed by O, V =1} IV T*/3 is the pressure, 
9a = -T* u ±x u " is the energy flow as seen by O and r Ap = T* u (-L^-Lp -g -Up-L^) is the 
spatial and traceless viscous-stress tensor. 5 The observed values of the fluid quantities are 
not independent of the choice of observer and are not just gauge effects. We discuss this in 
Appendix B. As implied by the decomposition of the second derivative of a scalar eq. (2.11), 
EMTs containing second derivatives will in general have contributions to all the components 
of the EMT, implying that the EMT has imperfect-fluid form: it is impossible to find a 
velocity frame in which it could be decomposed and possess a perfect-fluid form. 

At least for some models (in particular, the one we discuss as a worked example in 
section 3), we can rewrite the energy flow q\ and viscous stress t^ u in terms of a scalar 
potential with an appropriate number of derivatives: 

q x = ^ Vpg, (2.15) 

TAp = (l$± v p -^Ap^) V^V.tt 

with q the potential for energy flow and it the anisotropic stress. 

We present a detailed discussion of the conservation equations for the EMT using the 
language of covariant decomposition in Appendix A. Here we are just going to discuss the lin- 
ear scalar perturbations on a spatially flat Friedmann-Lemaitre-Robertson- Walker ("FLRW") 
metric working in the Newtonian gauge, assuming that vector and tensor perturbations be 
negligible, 

ds 2 = -(1 + 2^)dt 2 + a 2 (t)(l + 2<Z>)dx 2 . (2.16) 

The exact equations presented in the appendix, (A. 8) and (A. 9), can be perturbed directly. 
Doing so, we obtain the standard results. At the background level only the timelike equation 
contributes, giving the energy-conservation equation, 

£ + 3H(£ + V)=0, (2.17) 

5 The quantity qx is usually called the heat flow. However, as was demonstrated in [58] for kinetic gravity 
braiding theories, entropy production and heat flow are zero while qx does not vanish. In that case, the vector 
qx quantifies the energy carried by charges moving in the scalar frame. 
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with the overdot representing the derivative with respect to coordinate time t. At the linear 
level in perturbations, we can transform into momentum space obtaining: 



SS + :! /./ i AS + <vP) + I £ + V){<-> + <I> ) - ^ + 9© = 0, (2.18) 



/ k 2 \ k 2 

(£ + v)\e + 2He--~y\ --z6v + ev+ (2.19) 



a 2 / a 2 



where the 5 prefix signifies a perturbed quantity and we have defined the divergence of the 
perturbation of spatial velocity for the observer in the scalar frame: 6 

e = ikM = — % . (2.20) 

ma z 

It is customary to make the equations (2.18) and (2.19) dimensionless by dividing through 
by the background energy density: for us, it will prove more convenient to keep them in this 
form, since we do not wish to at this stage follow the usual practice of defining a sound speed, 
etc. and we wish to avoid any potential singularities this procedure might introduce in the 
variables. 



The form of the equations (2.18) and (2.19) is in fact gauge invariant, since in the 
Newtonian gauge eq. (2.16) the metric potentials correspond to the gauge invariant Bardeen 
potentials [15]. Gauge invariance, however, does not imply frame invariance. The standard 
references (e.g. ref. [16]) perform the decomposition of the EMT eq. (2.14) in the rest frame 
of the fluid and therefore do not contain the energy flow terms involving the potential q. 
The difference between any two frames is a Lorentz boost between the two velocities of the 
observers in whose rest frame the EMT is decomposed. Given that the background is FLRW, 
there is only one reasonable frame on the background level, that of the background comoving 
observer. Thus all the frames differ purely on the level of perturbations, with the boosts 
transforming equations (2.18) and (2.19) by trading the divergence of energy flow k 2 5q for a 
redefinition of the divergence of the momentum perturbation (£ + V)@. 

For a model with second derivatives in the EMT, the choice of any of these related 
frames is perfectly fine apart from the rest frame. In vacuum configurations, £ + V = 
and the divergence of momentum perturbation {£ + V)@ vanishes. Therefore the rapidity 
of the boost required to compensate a divergence of energy flow k 2 5q tends to infinity as we 
approach the vacuum equation of state. This means that close to vacuum configurations, 
perturbation theory in terms of the underlying scalar-field variables itself breaks down if we 
require the frame to be a rest frame. 

We have provided more detail for this discussion in Appendix B. Suffice it to say that 
the presence of the energy flow potential q in the perturbation equations (2.18) and (2.19) will 
be the origin of much of the novel behaviour described in this paper. Scalar-field theories 
more general than k- essence always contain such a term in the EMT decomposed in the 
scalar frame (2.2). Close enough to vacuum configurations, the coefficients of the terms in 

6 Note that there is a relative factor of a between & defined above and the variable 8 used by Ma and 
Bertschinger [16, Eq. (23b)], so that 8mb = aQ. This is a result of our using coordinate instead of conformal 
time. 
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(2.18) and (2.19) driving the evolution of perturbations under normal circumstances vanish 
and that role is taken over by the q terms. This then allows for such behaviour as a stable 
violation of the null-energy condition [47, 62] or "phantom crossing". The existence of the 
additional derivatives on the 5q perturbations compared to the G terms will mean that the 
transition between the two behaviours will not only depend on the equation of state of the 
dark energy but also on scale. 

Another aspect worth noting here is that the background value of q appears in eqs (2.18) 
and (2.19). This does not have to vanish, only its spatial gradient must at the background 
level (the total value of q is some function of the scalar field, <fi and m, so its background 
value is going to be that function of the background values of 4> and m). 

The aim of this discussion is to reduce the perturbation equations (2.18) and (2.19) 
into a single evolution equation for the density perturbation 5£. In the rest frame, this is 
simple since one can just eliminate the variable 0. In the scalar frame, in general both 
and 5q have to be eliminated. Luckily, one can exploit the structure of the conservation 
equations arising from their transformations under Lorentz boosts of u M to make an algebraic 
redefinition of to define the divergence of the total energy flow, 

3= (£+V)@-^-+q@. (2.21) 

3 is indeed proportional to the velocity divergence in the rest frame, ©ll as defined in 
eq. (B.6). However, since it is only an algebraic redefinition for us, the meaning of the 
perturbation variables remains: S£, 5V, etc. are still those measured in the scalar frame and 
not in the rest frame. Importantly, 3 is not singular when w = —1. This definition allows us 
to transform perturbation equations (2.18) and (2.19) into a simpler set, 

6£ + 3H(6£ + SV) + S+(£ + V)$> = 0, (2.22) 

~+bH~ + {£ + V)—„ 5V- 5tt = 0. (2.23) 

cr a z 3 a 4 

The above then need to be furnished with Einstein's equations. We obtain the independent 
equations by projecting out particular parts of the Einstein tensor and then perturbing, 

SiGauU^vT) =2^- + QH® - QHH = S£ , (2.24) 

2^ =-S, (2.25) 

$ + V = 5ir, (2.26) 

where eq. (2.25) makes the meaning of 3 as the total energy flow explicit. Contributions 
from external matter would need to be added to both the sets of equations above. 

In order to solve for the evolution of the energy-density perturbations, we need two 
extra ingredients: the closure relations between 8V, Sir and 8£. We shall first discuss these 
relations in perfect-fluid hydrodynamics and compare them with k-essence theories. We show 
that even for a canonical scalar theory with a potential the standard relation between 8V 
and 8£ is not hydrodynamical, despite the EMT having perfect-fluid form. This departure 
from perfect-fluid hydrodynamics becomes only worse for more-general scalar-field models. 
Since all EMTs obey equations (2.22) and (2.23) at linear level, the closure relations provide 
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the totality of information that can be accessed from gravitational interactions at the linear 
level. 

On the other hand, since we are investigating theories which are described by classical 
equations of motion, it is possible to obtain the exact closure relations. As we are dealing 
with equations of motion for a scalar field of second order, we can always express the closure 
relations in terms of variables defined on a constant-time slice. In section 2.5, we will present 
a general parameterisation which will serve us as a book-keeping device to ease these calcula- 
tions for particular classes of models. In section 3 we will take a specific model and perform 
this calculation explicitly. Our aim will be to obtain closure relations that to the largest 
extent possible depend on a set of physical parameters of the theory, such as the equation of 
state or the speed of sound. 

2.3 Perfect fluids and hydrodynamics 

In this section we will review the properties of perfect fluids. For a more detailed discussion, 
we recommend the review [75] . The perfect- fluid EMT is obtained by averaging a distribution 
of particles to obtain fluid elements and takes the form 

= S UfJL u v +?V, (2.27) 

with the energy density and pressure the only fluid variables. In Newtonian hydrodynam- 
ics one starts from the axioms of energy and momentum conservation, with the separate 
requirement that mass be conserved. From these the conservation of the EMT is obtained. 
This has to be appropriately generalised for relativistic fluids. The energy and momentum 
conservation equations become the continuity and Euler equations, while mass is no longer 
conserved. This axiom needs to be generalised to some appropriate charge conservation. 

The velocity field defines the rest frame (or the lagrangian observer) for the fluid 
elements. This velocity must be by definition time-like since in the construction of the EMT 
it is the tangent to the world-line of individual fluid elements. It is of course possible to 
observe this EMT in a different frame. However, it will then not in general have a perfect- 
fluid form and the measured energy density and pressure will be appropriately transformed. 
We discuss the issue of frame choice in detail in appendix B. In the rest frame, we can relate 
the pressure to the energy density through 

V = V(£,S), (2.28) 

where S is the entropy. In cosmology one often considers either adiabatic or entropy perturb- 
ations. The former are such perturbations of the configuration of the fluid that the pressure 
and energy-density perturbations are proportional to each other. The latter are such that 
despite a non-vanishing pressure perturbation, the energy density is not affected, i.e. these 
are perturbations only of composition of the fluid assuming constant energy density. The 
perturbations can be written as 

5V = c 2 J8 + a5S (2.29) 

where we have defined the speed of sound Cg as the coefficient in this equation. This sound 
speed then determines such physical properties as the Jeans length. As we have shown 
in appendix B, linear frame transformations affect the fluid variables at second order and 
therefore this form is independent of the choice of frame perturbation. 

For the purpose of calculating the anisotropies in the cosmic microwave background 
("CMB"), when the photons and baryons are strongly coupled and in equilibrium, one can 
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think of the photon-baryon mixture as making up one fluid. The entropy perturbations have 
a clear meaning of being a perturbation in the relative densities of photons and baryons while 
keeping the energy density constant. In such a case, the expression (2.29) is easily calculable 
[76, pg. 307]. This language is useful since a proof exists that at very large scales, the entropy 
perturbations are suppressed and only the adiabatic ones contribute [59, 77, 78]. However, in 
late-time cosmology, we are performing observations inside the Hubble horizon and therefore 
the entropy perturbations cannot in general be neglected. 

As we will show here, this description is not so helpful for classical scalar-field theories 
when applied to the late universe. The issue boils down to the question of what the degrees 
of freedom that are being perturbed in order to obtain the expression (2.29) are. There exists 
another description of fluids, much closer to our classical-scalar-field case, in terms of the 
thermodynamical potentials, temperature T and chemical potential [i. For a hydrodynamical 
perfect fluid, it is always possible to rewrite the system (2.28) in terms of these two potentials, 

V = V(T, f j 1 ), (2.30) 
£ = £(T, f x). 

The pressure and energy-density perturbations are now completely determined by the per- 
turbations in the values of T and [x. In general, we now have two types of propagating modes, 
one corresponding to temperature perturbations, the other to chemical-potential perturba- 
tions, i.e. effectively composition. 

Let us for the moment simplify and assume that the fluid be adiabatic: such a fluid is 
dependent on just one thermodynamical potential, e.g. the temperature, V(T), £(T). For 
the avoidance of ambiguity, let us also specify that we are observing this fluid in the rest 
frame. The perturbations of energy density and pressure in the rest frame are then related 
by 

5V = c 2 S£ , (2.31) 

where c 2 = V t T/£,T must have the same meaning as the sound speed defined in (2.29). In 
this language, the temperature is the degree of freedom which determines everything about 
the configuration of the fluid. Indeed, when discussing the photon gas, one does tend to talk 
of temperature perturbations rather than energy-density perturbations. For a photon fluid 
the chemical potential vanishes, since photons carry no charge, and thus everything depends 
just on the temperature. Perturbing the temperature perturbs the Bose-Einstein occupation 
function and therefore the number density of photons, the energy density and pressure. From 
this discussion, it should be clear that the evolution of the total pressure and energy density 
with respect to a time coordinate in the frame of the observer r is also just related by the 
same sound speed, 

V = c 2 J , (2.32) 

since everything is just a function of T and that c 2 = w, the equation-of-state parameter. 

In the non-adiabatic case, when both the temperature and the chemical potential are 
relevant, there exists another direction in which the system could fluctuate. The pressure 
perturbation can now be written down as 

SV = P iT ST + V tlt Sfi , 
5£ = £ t T^T + £^8/i . 
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We can also write down a very similar expression for the time evolution of the total quantities 



V = v, t t + P> , 
£ = £ <T f + £^fi . 

We can eliminate the 5T perturbation by combining these two sets of equations to obtain 

5V = C 2 5£ + (V- C 2 £) (2.33) 

with C 2 = V^t/£,t- We can now compare this expression to eq. (2.29) and see that they 
have the same form. We can read off the definition of the sound speed in terms of the 
thermodynamic potentials and find an expression for the entropy perturbation, in terms of 
5fj,. Note that we could have eliminated the chemical potential instead of the temperature 
perturbation. It is important to stress that the form of the expression (2.33) is independent of 
the frame of the observer, provided that it differs from the rest frame only at linear order. In 
particular it does not rely on coordinates, but is purely determined by the fluid configuration. 

The idea of this rewriting of the system (2.28) in terms of the thermodynamical poten- 
tials is that we can rewrite the system in terms of the degrees of freedom that we are free 
to perturb independently on any spatial hypersurface and then obtain an effective equation 
(2.33) eliminating some of them. As we will see, in a classical scalar-field theory, we also have 
such underlying degrees of freedom (the field value and the magnitude of its first derivative) 
on which all the fluid variables depend. The EMTs of the scalar systems look more similar 
to eq. (2.30) than to eq. (2.28). 

A fundamental assumption that we have used in the above is that no hydrodynamic 
quantity is an explicit function of time: pressure and energy both depend only on the ther- 
modynamical potentials. The potentials evolve according to the laws of thermodynamics, 
changing the configuration of the fluid, but if the potentials were to be returned to the same 
values as at some point in the past, the pressure and energy density would also return to 
their previous values. This is why the phase diagram of water does not change over time. 

One could imagine violating this principle. This would really mean that what the fluid 
is changes as a function of time, not just that its configuration evolves. Let us go back to 
a single thermodynamic potential, T, but this time together with an explicit dependence on 
rest-frame time r. We can write the derivative w.r.t. the rest-frame time r, 7 

f T=Br V + V,f T , (2.34) 
dr dr 

At a single point in time, i.e. on some spatial hypersurface in the rest frame, the perturbations 
will look as they do for the adiabatic fluid, since the only perturbed variable will be T, 

5V T f = V, T ST, (2.35) 
S£ r { = £ t&T . 



7 Note that using any other time coordinate here would mix the meaning of what should be considered to 
be a perturbation of the fluid in the rest frame and a perturbation of the rest frame itself. We should always 
rewrite all the variables to be in terms of rest-frame quantities where the fluid has a perfect form. 
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However now, when we transform into the frame of a different observer, the value of the 
perturbation observed by them is going to be corrected by the fact that the rest-frame time 
and the proper-time of the two observers are no longer the same, 

5V = d T V5T + 5V Ti = V, T (ST - T<5r) + V5t , (2.36) 
SS = d T £5r + 6£r{ = £,t (ST - TSrj + £5t , 

where we have eliminated the partial derivatives w.r.t. r using eqs (2.34). Using the definition 
of the sound speed in the rest frame, as we had defined it for the adiabatic fluid, we can 
combine the above two results to obtain, 

6V = c 2 s 5£ +{V- c 2 £^j 5t . (2.37) 

This relation has a very similar form to that for the non-adiabatic fluid, eq. (2.33), which was 
defined as being dependent on at least two thermodynamic potentials and with no explicit 
dependence on time. However, its meaning is completely different: here the substance depends 
only on one thermodynamic potential, but the way it depends on it evolves as a function of the 
rest-frame time r. Thus it is not the same fluid at two different times. Usual hydrodynamics 
does not apply in these circumstances. Or more bluntly put: the phase diagram of such 
substances would not be the same at different points in time. 

If we now assume that at the background level the rest frame and the observer's frame 
are the same, the transformation between the two frames is small. The relative perturbations 
of coordinates are linear. Since r is the proper time for the rest frame, we can think of the 
velocity of the rest frame as the derivative 

U„ = -d^r , (2.38) 

with the normalisation automatically taken care of. We can thus arrive at the standard 
expression [66] for the pressure perturbation of this substance 

5V = c 2 s 5£ + £(cl-c 2 )^, (2.39) 

where we have defined the "adiabatic sound speed" c 2 as V j£ and V = di5U l is the diver- 
gence of the velocity field perturbation. We should underline once more: this expression 
is not hydrodynamical. It explicitly depends on the coordinates of the frame in which the 
observations are being performed. V should not be thought of as an entropy perturbation. 

2.4 k- essence closure relations 

We have spent quite some time laying out the hydrodynamical picture in order to facilitate a 
comparison with a classical scalar-field theory. Here we will discuss general k- essence, defined 
by the Lagrangian 

C = K{<j>,X). (2.40) 
In the scalar-field variables, the EMT can be written as 

V = ( mK m ~ K) -^- 2 — + K [g^ + -^— 2 — ) , (2.41) 
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where we are using the definition m = y2X. We've written the EMT in such a way that 
allows us to easily see that it takes a perfect-fluid form when we identify the gradient of the 
scalar with the rest-frame velocity, i.e. = V«0/m. As we have already stated in section 2, 
this choice identifies <p with a comoving clock for the scalar system and therefore to = u^V^cj) 
with the rate of flow of this time. There is no other choice of frame we could have made for 
this model, since only with this choice of is the EMT of perfect-fluid form. The EMTs of 
more general scalar-field models cannot be rewritten in perfect-fluid form at all [58]. 

It is important to stress that this perfect-fluid description does not always apply for 
the EMT (2.41). As we have stated in section 2.3 the velocity needs to be time-like so 
that it can represent the velocity of fluid elements. A configuration of a scalar field can in 
general have non-time- like gradients (for example, a static domain wall). There is nothing 
wrong with such a configuration: its evolution is described by the equation of motion for 
the scalar field in the usual way. However, the EMT for such can not be interpreted as 
one of a perfect fluid. Similarly, the EMT of an oscillating scalar configuration is not of 
perfect-fluid form at all times. The velocity (j) vanishes and reverses direction. This means 
that at the rate of flow of time vanishes and reverses and the scalar reaches the maximal 
amplitude. In between two maxima, we can interpret the scalar as a perfect fluid, but this 
interpretation breaks down at the maxima themselves. Of course, it is well known that an 
oscillating scalar field (for example, at the end of inflation) effectively evolves as pressureless 
dust. Indeed this is true. However, there we are making a statement about the properties 
of the EMT averaged over oscillation cycles. This requires the existence of an external clock 
which describes the oscillatory solutions and therefore the existence of an external reference 
system which increases monotonically. The EMT of the scalar in such a situation is not that 
of a perfect fluid at all times, it is just that the particular solution can be described by an 
averaged EMT which has the form of the EMT for dust. 

Perfect-fluid evolution can always be written in terms of a non-canonical scalar-field 
theory [79]; however, not all solutions of k- essence theories can be described by perfect flu- 
ids. The reason why fluid descriptions are at all useful in cosmology is that the symmetry of 
the FLRW background requires that the spatial gradients vanish. Thus particular solutions 
always have time-like velocities (apart from the oscillatory cases mentioned). Linear per- 
turbation on this background by definition do not change this property since they are small. 
However, the question of whether the particular solution behaves like a fluid is in principle 
very different from the question of whether the scalar-field itself is a fluid. As we show here, 
this is not always the case and therefore it should not be a surprise that one does not in 
general obtain hydrodynamical relations between the perturbations of pressure and energy 
density. We will show this explicitly now. 

For ^-essence models, both the pressure and the energy density are functions of just 
two variables, (j) and to. As we discussed on page 13, it is natural to first write down the 
perturbations in terms of these underlying independent variables, 

sv = r m Sm + r^5<p , 

5£ = £ m Sm + £<p54> ■ 

where the subscripts m, 4> denote partial differentiation w.r.t. those variables. We can then 
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eliminate the derivatives w.r.t. <p by taking time derivatives of £ and P, obtaining 
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Given the definition of the spatial velocity divergence O, eq. (2.20), we can combine the 
above to obtain the k- essence closure relation 



where C 2 = P m /£ m . This relationship is not the hydro dynamical closure relation for a 
non-adiabatic fluid (2.33), but rather it is of the type for a substance depending explicitly 
on its rest-frame time (2.39), exactly since it depends on the perturbation of the rest-frame 
represented by O. Even assuming that the gradients V M are time-like, a general k-essence 
theory is not hydro dynamical, but represents a fluid which itself changes with time. There 
is only one thermodynamical potential, m, and there are no entropy perturbations in this 
model. However, in the shift-symmetric case, K = K(X), there is no explicit </> dependence, 
as a result of which c 2 = C 2 . In this case, the coordinate dependence of the closure relation 
(2.43) disappears and we recover the standard behaviour for a perfect adiabatic fluid, eq. 
(2.31). However, this is a fluid which can only carry scalar perturbations and not vorticity: 
a superfluid. 8 

One final aspect that needs to be checked is the relation of C 2 to what we call the 
physical sound speed c 2 . C 2 relates the pressure and energy-density perturbation through 
eq. (2.43), and therefore determines the Jeans length. What we will call the physical sound 
speed Cg in the case of a scalar-field theory is defined as the speed of propagation of small 
perturbations of the field. This physical sound speed determines the rate at which information 
propagates and therefore the causal structure of the scalar medium [80]. It is obtained by 
calculating the effective metric for perturbations. We show how this can be done in section 
3.2. In general the two quantities do not have to be equal, but, in the case of k-essence 
C 2 = for both shift-symmetric and general models. The Jeans length in k-essence is 
therefore always determined by this physical sound speed c 2 , rather than any other speed 
of pressure waves [81]. This also could have been a different result and indeed it is highly 
non-trivial to obtain it in more complicated theories. However, because of the link between 
the Jeans term and causality, resulting from the fact that it describes the reaction of the 
fluid to changes in the configuration, it is this physical sound speed of field perturbations 
that will in the end determine the behaviour of the scalar fluid at small enough scales. 

In conclusion: k-essence models hydrodynamics of an adiabatic fluid only in the shift- 
symmetric case. Away from shift symmetry, the substance it models changes its properties 
as a function of time. However, it still depends only on one thermodynamical potential 
and does not carry entropy perturbations. One can see that a non-adiabatic perfect fluid 
would be modelled by two shift-symmetric scalar fields. The two magnitudes of the gradients 
would represent two different thermodynamic potentials, while the independence from the 

8 As we have already stated, this is a result of the Frobenius theorem. The rotation tensor for the velocity 
vector field vanishes, u; M „ = V[ M u„j = 0, since u M is a derivative of a scalar. The vorticity vector for the 
vector field u M is then defined and therefore also is zero. 




(2.43) 
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field values would ensure that the fluid were unchanging. The description of perfect-fluid 
hydrodynamics through scalar theories was discussed in ref. [79] . 

As we will start to show in the next section, the closure relations for more complicated 
models depart even further from the hydrodynamical picture. Nonetheless, just as in the 
case of eq. (2.42), hydrodynamic quantities such as the equation of state and the sound 
speed do appear in the closure relations and therefore one can hope to express the general 
behaviour using a small set of physical parameters. We shall therefore be very careful when 
deriving the closure relations and will show explicitly that the correct relations are not those 
of hydrodynamics. 

Let us invest a little time to rewrite the perturbations for k-essence (2.42) using a 
different set of variables, which will make the transition to more difficult cases discussed in 
section 3 more transparent. The pressure perturbation can be considered to be a sum of two 
parts: spatial and temporal, 

k 2 

—5V = mV m A + V® . (2.44) 
ar 

with a similar result for the energy-density perturbations. We have defined the variable A, 

A=(^-™e)=e + 2He-^, (2.45) 

\a z mml a A 

to represent the perturbation of m as it would be seen in the scalar frame u^, i.e. it is 
the equivalent of the rest-frame perturbation (2.35). It is also equivalent to the spatial 
divergence of the perturbation of the acceleration a\ (see eq. (2.7)). Spatial in this case 
implies independent of the perturbation of the clock <p. 

We then obtain the closure relations by rewriting the velocity divergence in eq. (2.44) 
through the total energy flux 3 (eq. (2.21)) and then eliminating A from the pressure 
perturbation equation using the energy-density perturbation. We have to do this since only 
S and not A has an evolution equation (eq. (2.23)) in the language that we have set up. In 
this language the closure relations (2.43) for k- essence become 

k- essence: SV = c 2 5£ + 3(<£ ~ c l) (^^j J » ( 2 - 46 ) 
k 2 

—8tt = 0. (2.47) 
a z 

The horizon suppression factor, aH/k, appears here since S was defined in eq. (2.21) as 
containing two spatial derivatives already and this factor cancels out these derivatives. The 
k- essence model contains just functions of <j> and m and does not contain any higher-order 
derivative terms, therefore the perturbations themselves cannot contain such contributions. 
This changes in theories containing second derivatives of the scalar in the Lagrangian. Apart 
from this, the coefficients are just functions of time, dependent purely on the evolution of 
background quantities. It is only such terms that can reasonably be hoped to be constant or 
to scale simply when the scalar is following a scaling solution. Our aim in the next section 
is to extract this sort of behaviour for more general models. 

2.5 General parameterisation of closure relations 

As we have already mentioned, in order to properly solve for the evolution of energy-density 
perturbations in any particular model in the fluid language, we need to calculate the ap- 
propriate closure relations for it. At this stage it is useful to introduce a parameterisation 



-17- 



for these relations. This parameterisation will be based on the expectations of the kind of 
terms that more general theories of a single scalar degree of freedom yield. Our aim is not to 
introduce a set of parameters that will span the space of all possible and impossible theories 
and which should be constrained from the data. After all, the most general theory of a single 
scalar degree of freedom is already known [49, 50]. What remains is to perform calculations 
for these healthy theories within some a physically meaningful framework. Our parameterisa- 
tion should be thought of as a book-keeping device which should aid in these calculations. 
As we will show, the advantage of this approach, as opposed to just solving the equations 
of motion for the scalar fields, is that it allows for understanding the effect of changing the 
hydro dynamical properties such as the sound speed and equation of state, within particular 
classes of models, without worrying about the exact form of the Lagrangian. 

We propose that, under the assumptions we detail below, the form of the closure rela- 
tions will be contained within the structure: 

SV = CHE + 3Ei 2 § + S 2 | + foe-*6p + 7 P^e m (^j * , (2.48) 

-57T = U5£ + 3w 1 ^—j - + ro2 - + /V ^ + 7^e m ^-j . 

where p and 5p are the energy density and its perturbation of any matter external to the 
dark energy assumed here to behave as dust, O m is the velocity divergence of the external 
matter, e K is the effective Planck mass which appears as a result of our having eliminated 
Guv from the EMT of the dark energy. The set of parameters which need to be calculated 
for a particular class of models are: 1,2, ^1,2, Pv,-k an d 7P,7T- 

Let us justify this structure. Firstly, it is clearly not a structure which would be given 
by hydrodynamics of a perfect fluid, just as it wasn't in the case of ^-essence. In general, the 
EMT will contain standard k- essence contributions, which will result in the type of terms 
obtained in eq. (2.46). More general scalar models contain second derivatives in the EMT, 
generating entries in all the components, as can be seen from the decomposition (2.11). The 
fluid variables can for example depend on the expansion, 6 or the intrinsic curvature of the 
spatial hypersurface, ^>R. When perturbed these give 

69 = 6 + 3($ - HV) , (2.49) 

(5 (3) jR = 4 ^_ 4i?e . 

a 1 

The gravitational potentials here can be solved for by using the Einstein equations (2.24), at 
the cost of introducing a dependence on the energy density and velocity divergence of any 
external matter present in the cosmology. Thus in general models with second derivatives, 
the perturbed fluid variables will not only be a function of the scalar configuration variables 
84> and 5m, but also be functions of the configuration of the external matter. This is the 
origin of the maybe surprising parameters /3-p i7r and 7p i7r in the parameterisation (2.48). This 
doesn't make the system unsolvable, since there will also be conservation equations for the 
external matter equivalent to eqs (2.22) and (2.23), but it makes these more-general models 
significantly more complex, eqs (2.49) show that the dependence on O arising from the 
second-derivative terms is unsuppressed by k/aH, leading to the contributions of terms E 2 
and W2 to the parameterisation (2.48). 

Another general feature that appears in these models with second derivatives is the 
presence of dm /dr, the second time derivative of the scalar field, in quantities such as V (see 
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e.g. eq. (3.8) or ref. [58]). This is yet again another independent variable on which the 
pressure perturbation will depend. The method we propose for solving this complication is 
to use the equation of motion of the scalar field to eliminate dm /dr, prior to perturbing. Since 
any coupling to external matter will by definition also appear in this equation of motion, this 
sort of replacement will again give a contribution to the pressure of dark energy EMT from 
the external matter. In the equation of motion this coupling can be more general than just 
to the energy density of the external matter. 9 We will neglect the pressure of the external 
matter in this paper and thus have not included it in eqs. (2.48). 

In section 3 we illustrate the above somewhat abstract arguments with a concrete ex- 
ample model. Even this simplest example will require most of the structure of the para- 
meterisation (2.48). We will see that theories with second-derivatives in the EMT feature 
a new scale for perturbations: it is determined by the relative sizes of the k-essence and 
second-derivative terms in the perturbed EMT. At large enough scales, the k-essence terms 
will dominate and the perturbations will evolve as they do in k-essence. At smaller scales, 
the two-derivative terms will dominate bringing new features to the evolution of the per- 
turbations, such as anisotropic stress. This means that the parameters introduced in (2.48) 
are not strictly speaking functions of just the background. They will contain exactly this 
perfect /imperfect transition scale, and will in general interpolate between two values either 
side of this scale. However, the way that we have defined them, away from this transition 
scale they will behave as constants in k (not necessarily in time). In particular, no higher 
derivative terms are possible (higher powers of k /aH) if we limit ourselves to models with no 
more than two derivatives. We should also note that each class of Lagrangians with second 
derivatives will introduce its own transitions scale, dependent on the value of the coefficient 
it carries. In this paper, we shall only investigate a toy model with one such term. 

We should remark that this transition between the two limits is in the end determined 
by the relative importance of the and 5q terms in 3, eq. (2.21): the dominance of the 
Sq terms signifies the dominance of the two-derivative terms. It will clearly occur at small- 
enough scales, when the mode k is high. However, as the form of S, eq. (2.21), shows, the 
q terms will also dominate whenever the equation of state is close enough to that of the 
vacuum, £ + V = 0. Thus in models which would fit the observed expansion history, having 
w ~ — 1, the transition to the imperfect behaviour will lie at relatively large scales even when 
the two-derivative corrections are relatively suppressed and do not contribute significantly 
to the evolution of the background. It is interesting to note that this transition scale is not 
present in the equation of motion for the scalar. It is purely a result of the structure of 
the EMT. However, since it is the EMT that determines the behaviour of the gravitational 
potentials, this scale is capable of having an observable effect. 

We can now insert the parameterisation (2.48) into the perturbation equations (2.22) 

9 e.g. in Brans-Dicke theories, the equation of motion couples the scalar field to the trace of the external 
EMT, i.e. p — 3p, but it is a different combination of pressure and energy for galileon terms. 
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and (2.23) to obtain 



S£ + 3#(1 + C 2 )5£ +( 1 + 3S 2 + 9E X ( — ) ]S= (2.50) 
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ff* ">-{3 S 2 + r2 ]]ff=- (2.5!) 
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where we have explicitly included the non-conservation terms Sip which appear as a result 
of the evolving Planck mass in the presence of external matter, eq. (2.13), 

S 1 = 5 U-*T%*u>'V v >e\ , (2.52) 



where the 5 signifies that the term inside the parentheses must be perturbed. We have 
also dropped some terms involving the gravitational potential, subdominant subhorizon (i.e. 
negligible on scales k 3> aH, which are the most interesting ones from the observational 
point of view). Given this form, we can combine the two equations to eliminate the variable 
H and obtain a single second-order differential equation for 5£ with sources related to the 
presence of external matter. Let us make some simplifications first: 

• The term in eq. (2.51) containing £2 + \ w 2 wm contribute to the friction term in the 
evolution equation for 5£. When rewritten in In a as the time coordinate, it can be 
easily seen that this will be a contribution exponentially growing in time, large at small 
scales. Depending on the sign, it either is an anti-friction term leading to small-scale 
instability or a very quickly growing friction term, destroying any perturbation growth. 
In the model we describe in (3), it is never important and vanishes on small scales. We 
suspect that this is in general the case for all consistent theories. We will set it to zero 
in what follows. 

• Since the pressure and anisotropy perturbations contain the external matter perturb- 
ation, Sp, the Jeans term in eq. (2.23) gives a scale-dependent source for the DE per- 
turbation. We can estimate from eq. (2.51) that, at small enough scales, the particular 
solution of the differential equation for 5£ we will approximately have is 

S£ - when + l^C* + 2 -H » (^) 2 . (2.53) 

These will be the solution provided that the homogeneous modes decay quickly enough. 
This is usually referred to as the quasi-static limit. Since in this configuration the 
dark-energy energy density tracks the dark matter one, this would be interpreted as a 
modification of the effective Planck constant in the Poisson equation, or a fifth force, and 
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is the sort of solution taken as the f(R) modified-gravity regime. The rather interesting 
observation is that in models with a DE/DM coupling, it is on scales smaller than the 
Jeans length that we observe such effective modifications of the Newton's constant. It is 
below the Jeans length that the dark energy is clustered in models where it is coupled to 
the DM. However, we must be careful as to whether this quasi-static solution neglecting 
the solutions to the homogeneous equations is in fact dominant. This may not always 
be the case. 10 

Neglecting the external energy density for the moment, we can obtain an evolution equation 
for the energy-density perturbation of the dark energy when it is isolated: 

+ 8 + 3((7 2 -£i- jW) + ^w 2 )s£+ (2.54) 



+ (1 + 3£ 2 ) ( C z + -II J -2 5£ + H 



iL( 4 (i_ii)-3£ 2 (i-6n)) + 



+ (3 - 2IT) 1 5 + ^w 2 - 3 (Si + ^-wx 



5£ = 



where we have dropped terms suppressed below the horizon and assumed £ 2 and C 2 vary 
slowly. 

At this stage let us stress that it is the combination of parameters (1 + 3E 2 )(C 2 + 2 /3H) 
that determines the Jeans length. As we will see, in our example model, it is this combination 
that will reduce to the physical sound speed squared, c 2 at small enough scales. Thus both the 
Jeans length and the speed of propagation of perturbations of the scalar will be determined 
by the same physical sound speed rather than a simple parameter relating 5V and 5£. 

We would also like to draw the reader's attention to the coefficient of the friction term 
in eq. (2.54). It contains corrections tui j2 which appear whenever the anisotropic stress 
is present. When the fluid is imperfect, the friction coefficient is changed (in our example 
it is reduced) and this modifies the power law with which the homogeneous modes of eq. 
(2.54) evolve, making 5£ decay more slowly (and the density contrast grow significantly 
faster) compared to its behaviour in a perfect regime. In hydrodynamics, viscous stresses 
act to increase the friction, so this result is very much against that intuition, underlining the 
non-hydrodynamical behaviour of these scalar models. 

We will now turn to a discussion of a particular model that features second derivatives 
in the EMT within the framework proposed in this general discussion. 



3 Worked example: non-minimally coupled k- essence 

The (re-)discovery of galileon-type theories [43-50] has significantly extended the model- 
building toolbox to include theories which describe imperfect fluids and can stably violate 
the null energy condition [58], allowing for the first time the study of models which have such 
phenomenology in a consistent manner. The fundamental difference as to more phenomen- 
ological approaches, is that these theories do have a well-defined Lagrangian and equations 

10 This is, at least in spirit, similar to the "adiabatic instability" in coupled models [82-85]. This instability 
is a result of assuming that the scalar field follows the minimum of an effective potential adiabatically, i.e. 
its evolution is quasi-static and is determined by the density perturbation of dark matter. However, there 
can exist circumstances when that is not the solution which the scalar obeys which is signified by a negative 
effective sound speed squared of the combined DM/DE fluid. 
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of motion which allow for a precise understanding of the evolution. Indeed cosmological 
perturbation theory in these models has already been extensively studied [86-97]. However, 
in all these models, the presence of second-order derivatives in the EMT complicates matters 
very significantly, obscuring the physical meaning of the results. 

The purpose of this paper is to show how one can use intuitive hydro dynamical lan- 
guage to study these more complicated theories. However, in the interest of simplicity and 
familiarity, we will study a much simpler model than the full Horndeski-type Lagrangian [50]: 
a k- essence theory non-minimally coupled to gravity. Again, this class of models has been 
studied before [98, 99], however from a very different point of view. We will focus on the fact 
that in such theories second derivatives in the EMT are generated through the non-minimal 
coupling to gravity. As a result of this, in the Jordan frame, this scalar theory can be seen to 
exhibit many of the properties of general galileon EMTs at the level of linear perturbations 
on a FLRW background. As we will see, this theory also features a rest frame different to the 
scalar frame as well as non-zero viscous stress thus addressing most of the possible structures 
in the more general theories. This worked example will allow us to demonstrate explicitly 
how one goes about deriving the parameters of the closure relations (2.48), resulting in a 
description which depends on physical properties of the fluid rather than a particular choice 
of Lagrangian to which the equation of motion is sensitive. 

As we have already mentioned in section 2.2, these models are usually described in the 
Einstein frame, where the coupling to gravity is removed at the price of a coupling between 
external matter and the scalar. In the Einstein frame, the second-derivative terms in the 
EMT are not present and the DE fluid is just a standard k-essence coupled to external 
matter. We choose to remain in the Jordan frame, since we are using this worked example 
to illustrate the features of more-complex models, such as kinetic gravity braiding, where no 
frame redefinition exists which would remove the second-derivative terms. 

3.1 General properties of energy-momentum tensor 

Our worked example is described by the following action: 



The first term defines the gravitational part of the action, which is the usual Einstein- 
Hilbert term but non-minimally coupled to a scalar <p through the function The 
second term, K(cf),X), is the usual k-essence term: a non-canonical kinetic term for the 
scalar field. Further, through £ ext we have indicated the possibility of some matter content 
external to the scalar field, which we do not consider in detail in this paper. Mp\ denotes 
the fundamental Planck mass, implying that the effective Newton's constant is 



We will use units in which Mpj = (87tGn) _1 = 1. The always-positive form of the non- 
minimal coupling, e x , ensures that the effective Newton's constant does not change sign. If 
we start the system off in a configuration where the gravitational sector does not have a 
ghost, it will never evolve to destabilise it. 

We will simplify the notation somewhat by using the overdot to signify the differentiation 
w.r.t. the scalar-frame time, / = d//dr in section 3.1 and 3.2. Once we go back to the 
discussion of perturbations in section 3.3, the overdot will again mean the differentiation 
w.r.t. background coordinate time in FLRW. 




(3.1) 



Cefr = G^e 



(3.2) 
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The Einstein equations are 

Gfa> = + e , (3-3) 

T* = e^K^V^V^ + e-*Kg^ + e"* (V M V, - ^D) e*W . (3.4) 

where is what we will refer to as the EMT of the dark energy and which includes the fact 
that we have divided the equations through by the effective Planck mass (see the discussion 
on page 7). T^f is EMT for any external matter. If we consider it, we will take it to be a 
perfect fluid, 

Tflu 1 = (Pext + Pext) Vfj,V v + Pextff/xi/ > (3-5) 

and we will assume that it be conserved, VT^ = 0. As was alluded to in the introduction 
to section 3, the non- minimal coupling to the Ricci scalar in the action (3.1) results in second 
derivatives appearing in the energy momentum tensor eq. (3.4). In this respect this model 
shares many features with the more complex theories featuring galileon terms or kinetic 
gravity braiding from the perspective of cosmological perturbation theory. The price of the 
varying Planck mass is the non-conservation of our choice of the DE EMT in the presence of 
external matter, 

= -xe-% c *V, (3.6) 
where the velocity is that defined in eq. (2.2). and 

k = mx^ , (3-7) 

where the subscript represents differentiation with respect to that variable. We can now 
decompose the dark-energy EMT in the scalar frame (2.2) according to (2.14): 

£ = e~ H E -kO, E = mK m - K (3.8) 

2f 2n • (™> 2 



T> = e -*p + m \x H + x*) + k ( — + -0 J , P = K 

qx = ka\ = V\q , q = -k 

Tpv = —ka^u = na^y , tt = —x 

The two variables E and P are the energy density and pressure arising directly from the 
k- essence term in the action, K((f>, X), with the subscript m signifying partial differentiation 
w.r.t. m. In the spirit of eqs. (2.15) we have explicitly provided the scalar potentials for the 
energy flow and viscous stress, q and tt. 

The decomposition (3.8) explicitly shows up the imperfect corrections to the EMT, all 
carrying two derivatives of the scalar <f>. In particular, we have a contribution to the energy 
density proportional to the expansion, a term in the pressure proportional to m and a non- 
vanishing energy flow in the scalar frame, proportional to the acceleration. These are the same 
modifications as in the case of kinetic gravity braiding [58], which provides the motivation 
for us to study this otherwise well-known model. In addition to the above corrections, in 
non-minimally coupled scalar models, non- vanishing viscous stress appears, here in the form 
of shear viscosity. There is also a bulk viscosity term in the pressure proportional to the 
expansion 6. 11 All these second-order terms appear with a common coefficient, k, which 



11 Since we are starting with a Lagrangian, the system must be conservative. It should therefore be stressed 
that these viscosity coefficients are not the same as one would obtain through the gradient expansion. They 
do imply that the evolution of the system depends on the shear and expansion, similarly to the way that shear 
and bulk viscosity do, but is not in fact dissipative. See the discussion in ref. [58, §3.7] where such a situation 
is discussed in the case of the kinetic-gravity-braiding model. 
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effectively provides a single additional parameter controlling such modifications. The model 
reduces to k-essence plus a constant correction to the Planck mass in the limit a = const. 

3.2 Equation of motion and effective metric for perturbations 

The equation of motion for the scalar that is obtained directly by varying the action (3.1) is 

V„ (K X V^) -K 4> = l -x^R . (3.9) 
where we obtain the Ricci tensor from the Einstein equations (3.3), 

R = S - W + e~ x (pext - 3p cx t) • (3.10) 

Combining these and rewriting the second derivatives of the scalar in terms of the kinematical 
decomposition (2.11) we obtain for the equation of motion, 



(E m + ^m^eA m+m (p m + ^m^e^ 9 + mE^ = (3.11) 

= (e - 3P - 3m 2 e z (x^ + + p cxt - 3p ex t) • 

Because of the presence of the second derivatives in £ and V and the non-minimal coupling, 
there is now a Planck-mass-suppressed contribution which corrects the coefficients of the 
second-derivative terms, m and 6. Using this form we can also explicitly see that this scalar 
is coupled to the trace of the external EMT. This has an important impact on the dynamics 
of the scalar in the presence of external matter. 

An important feature of the form of the scalar equations of motion (3.11), with all the 
second-derivative terms explicitly collected together, is that it allows us to calculate the sound 
speed by identifying the effective metric for the propagation of linear perturbations. The idea 
is that we linearise the equation of motion around some background configuration (including 
gravity) and then collect all the highest derivative terms. We will obtain an operator Quv, 

Q^VpVM + O (V„(ty, S4>, $, *) = . (3.12) 

Now given the standard eikonal ansatz, 50 = A(x) exp (iuS(x)), with some slowly varying 
amplitude A{x) and then taking the formal limit ui — > oo, we obtain 

a^V M <SV„<S = . (3.13) 

The tensor is now the effective contravariant metric for the propagation of perturbations 
[80, App. A]. In the case of our equation of motion (3.11) on an arbitrary background, it 
remains a simple diagonal metric, just as in k- essence, 

G»v = + u,u v + + 3 -,y) V . (3.14) 

Note that this is an exact expression valid for any background metric, not just, for example, 
for FLRW. We do have a correction to the k- essence result, with the non-minimal coupling 
to gravity resulting in Mpj-suppressed terms coming from the second-derivative terms in the 
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DE EMT. From this result, we can immediately read off the condition that the scalar degree 
of freedom not be a ghost, 12 

D = ^ + §x}e->0 > (3.15) 
m 2 Y 

which comes from the coefficient of the time-time part of the acoustic metric. In addition, 
the sound speed for the propagation of small scalar-field perturbations is the ratio of the 
diagonal space components of the metric (3.14) to the time-time part, i.e. 

4 = ^ + 3<3>0, (3.16) 
«*."> = 

where we have defined a new positive coupling f5 arising from the non-minimal coupling 
to gravity. Again, the sound speed squared needs to remain positive in order to prevent 
short-timescale gradient instabilities from developing, since this is the quantity that enters 
the dispersion relation. We can rewrite (3.16) in the better-known form, 

c2 = K X + |x 2 e" 

° s K x + 2XK XX + |x 2 e* ' 



2K X + 4XK XX + 3xje> 



The sound speed is a combination of the standard expression for k-essence and a correction 
from the non-minimal coupling. The sound speed is always equal to the speed of light for 
a canonical kinetic term. On ghost-condensate-type attractors (P m = 0) the sound speed 
is non-zero as a result of the modification from the coupling to gravity. Thus if the non- 
minimal coupling to gravity is weak (x(0) is slowly varying) and we have approximate shift 
symmetry in the K (0, X) term, we would expect the solution to approximately follow the 
k-essence attractor, with w ~ — 1 while the sound speed would naturally be close to 3/3, at 
least when the dark energy is dominant. 

In general, the sound speed c 2 > 3/3. This restriction could perhaps be avoided for some 
functions K((f>, X) involving inverse powers of X. However, these are likely to have additional 
pathologies in their phase space and we shall not consider such models here and assume this 
hierarchy. 

It should be stressed that the result (3.16) is valid for the action (3.1) and is not a 
general prescription for the sound speed. See e.g. refs [47, 100] for the modifications to 
the expressions generated by galileon-type terms in the Lagrangian. For the avoidance of 
confusion with other uses in the literature, we will refer to eq. (3.16) as the physical sound 
speed. The key remark here is that this metric (3.14) describes the effective space time (and 
its causal structure) that a perturbation in the scalar field sees. Thus signals sent using the 
scalar field would propagate in this metric and therefore with the speed (3.16). 



One should ask about the relevance of the speed of propagation of scalar signals to 
what is normally considered in the study of cosmological perturbations: the sound speed is 
normally taken to relate the pressure and energy-density perturbations; frequently it is even 



12 This is in addition to the condition G e ff > 0, which is required to keep the gravity sector healthy. 



defined as 5V/5S. One speaks about the speed of pressure waves with the idea that the 
response from pressure is what stops collapse and gives the Jeans length. As we will show, in 
the model (3.1) this reaction to collapse is not just driven by the pressure V alone, but also 
by the anisotropic stress, as was already previewed in eq. (2.54). Moreover, we will show that 
despite the fact that the relationship between the pressure and energy-density perturbations 
for scalar-field theories is in general complicated, the Jeans length is determined by the actual 
physical sound speed (3.16). Thus the Jeans length, even in these more complex theories is 
determined by causality and the rate of propagation of signals in the scalar fluid. There is only 
ever one physical sound speed and it is determined by the configuration of the cosmological 
background and therefore, in FLRW cosmology, can at most be a function of time. 

As we show in section 3.3, the quantity (3 is associated with the existence of a new scale 
delimiting two behaviours of cosmological perturbations: it separates perfect and imperfect 
regimes. We will show that this is a separate scale from the Jeans scale discussed above and 
does not appear in the equation of motion. It only arises as a result of the structure of the 
EMT. 

We shall now return to the discussion of perturbations in the fluid language. As we 
have stated in section 2.5, we need to eliminate any terms in the EMT which are not defined 
on the spatial hypersurface. In our model (3.1), this means eliminating rh from the pressure 
in the EMT (3.8). It is easy to rewrite the equation of motion (3.11) in a form which will be 
directly usable: 

= -kclO - ^ + Pe~ K (E-3(P + mV(x # + xj)) + p cxt - 3p cxt ) . (3.17) 

We can now use eq. (3.17) directly in eq. (3.8) to re-express the pressure as 

V =(1 - m (e"*P + m 2 (x# + 4)) + Pe~* (e - (3.18) 

- - jjj k8 + f3e~ K (p ext - 3p exp ) • 

We now see some of the structure which we previewed in section 2.5: the coupling of the 
scalar to external matter appears explicitly in the pressure with the coupling parameter f3. 
Secondly, a term proportional to the physical sound speed appears as a coefficient of the 
second-derivative term, 9. It is this term that will dominate the pressure at small scales and 
determine the Jeans scale. 

Another implication of the the pressure's containing m is the presence of singularities in 
the phase space: as can be explicitly seen in eq. (3.18). Whenever D — > 0, we will also have 
P — > oo and so a number of terms in the pressure will diverge. For a cosmological background 
solution, if the pressure in this limit is negative, the trajectory will evolve toward a pressure 
singularity: H will diverge within a finite time, while a and H remain finite [101]. This sort 
of behaviour can also be seen in phase spaces of models with kinetic gravity braiding in ref. 
[62] and arises in the same way. We do not observe this in minimally coupled k- essence, 
despite the possibility of having configurations with a vanishing D since the pressure does 
not contain m directly. 
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Let us close this discussion of the general properties of the models (3.1) by discussing 
the limit to f(R) gravity, which is a subclass. The action is 

S = J d A x^M^f(R) . (3.19) 

The standard procedure is to start off from this action and rewrite it in the Einstein frame, 
which explicitly shows that an extra scalar degree of freedom is present. However, even in 
the Jordan frame one can perform a Legendre transformation introducing a scalar explicitly 
and obtaining the action 

S = J d 4 xV^gM^ [U{4>)R + /(</>) - 4>U(4>)} . (3-20) 

equivalent to (3.19) [5]. This form of the action shows that f(R) gravity is equivalent to our 
class of models (3.1) in the limit where the k-essence function K is purely a potential V(4>). 
In the Jordan frame, the kinetic term for the scalar is purely generated through the coupling 
to gravity. In the EMT (3.8), apart from the potential energy V(4>), there are only terms 
arising from the second derivatives and the behaviour of the fluid is always determined solely 
by them. In particular, we have 

no scalar ghost D > 44> > , (3-21) 

c 2 s = 1, 
P = 1/3. 

Thus f(R) is a special case of the class of models (3.1) where the Jeans scale and the 
transition scale associated with /3 are unobservable, since Cg and /3 are both order one. In 
cosmology, both in the background evolution and in the perturbations, the second- derivative 
terms dominate the DE EMT on all scales. One can reverse this observation and note that we 
approximately recover the behaviour of f(R) solutions, including the properties (3.21) for any 
Lagrangian belonging to the whole class (3.1) whenever the background configuration is such 
that the kinetic energy arising from the k- essence term is small compared to the contribution 
to the energy density arising from the second-derivative terms, i.e. for background cosmology 
when 

e _x (E + P) < kH . (3.22) 

Since H contains a contribution from external matter, given everything else be fixed, this 
condition is more likely to be satisfied in the past, when dark matter dominates the universe, 
although whether this happens is solution dependent. This type of solution was studied 
under the name of the "strong-coupling" regime for Brans-Dicke theories in ref. [102]. 

When condition (3.22) is not satisfied implying that the k-essence terms dominate, the 
phenomenology of the theories (3.1) is much richer: both the sound speed and the coupling 
f3 can take small values, placing new potentially observable scales affecting the evolution of 
dark-energy density perturbations inside the cosmological horizon. In general, the history of 
the evolution of such a scalar will go through both the f(R)-like regime and this k- essence- 
dominated regime. The history of this evolution is highly model dependent. Since the purpose 
of this paper is to propose a prescription for studying such models, we will only illustrate 
our method in the somewhat unphysical case of a cosmology of an isolated dark-energy fluid, 
corresponding to the far future, when the dark matter has been diluted away by the expansion 
of the universe. This will allow us to demonstrate a number of general properties for scalar 
theories involving second derivatives in the EMT. However, we will return to a more detailed 
study of the full richness of these models in the presence of dark matter in a separate work. 
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3.3 Cosmological perturbations 

In this section, we will study the behaviour of the non-minimally coupled k- essence model 
(3.1) at the linear level in cosmological perturbations using the prescription introduced in 
section 2.5. For simplicity, we will only consider the case of a universe containing solely the 
scalar field. This will be sufficient to prove a number of statements we have made so far, 
but we will not explore the full richness of the models in this work. However, at appropriate 
moments, we will restore the dependence on the external matter to illustrate generic prop- 
erties. We will not address the issues of non-linear perturbations and the existence of the 
chameleon mechanism, which can change the behaviour of the scalar at small-enough scales. 

We will explicitly show that this model features a new length scale which divides the 
regimes of perfect and imperfect behaviour of the system. In the limit of f(R) theories, this 
length scale, as well as the Jeans scale, lie beyond observable scales. It arises in our fluid 
language as a result of the dominance in the pressure and energy-density perturbations of 
terms containing second derivatives of the scalar field. Such terms are generic in other models 
containing second derivatives, such as the galileon or kinetic-gravity braiding, therefore we 
expect such a length scale and transitions between two regimes with different perturbation 
behaviour to also be generic. Secondly, we will prove that the Jeans scale is purely dependent 
on the speed of propagation of the scalar-field perturbations, the physical sound speed (3.16) 
and not on the naive relationship between the pressure and energy-density perturbations. 

From here on now the overdot will signify the derivative w.r.t. the FLRW background 
coordinate time, d/dt. The prime represents the derivative w.r.t. In a. 

We will now explicitly demonstrate how the closure relations in our example model 
reduce to the form proposed in the parameterisation (2.48). Let us start by reminding 
ourselves of the expressions for the DE energy density and pressure, eqs. (3.8): 

£ =e-*E-^6, (3.23) 
dr 

V =e~» ((1 - 3/3) (P + m\x H + + P (e - ^ ) - (c? - |) 9^ . (3.24) 

These are exact covariant expressions in the absence of external matter and where the de- 
pendence on dm/dr has been eliminated through the equation of motion (3.17). In order to 
be explicit, let us define the equation-of-state parameter and the adiabatic sound speed, 

(3.25) 



3(1 + w) ' 

where the prime ' denotes a differentiation with respect to the logarithm of the scale factor, 
In a and all the quantities above are background quantities. 13 The w of eq. (3.25) is the 
equation of state for the DE; however, in the DM-less scenario being discussed here it is also 
the total w of the cosmological expansion. 

13 Note that the relationship between and w is strictly speaking only valid when no external matter is 
present: in the presence of external matter, the DE energy momentum tensor defined here is not conserved 
(see eq. (3.6)) and therefore £ has a new contribution dependent on p c *t- This modifies the expression for c^. 
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We will focus our study on configurations where 

< 1, (3.26) 

The first condition implies that the Planck mass be evolving slowly. The second implies that 
the external matter, which is universally coupled in this model, be coupled relatively weakly, 
so as to avoid strong violations of the equivalence principle. It also means that we are not 
in the well-studied /(i?)-like regime discussed on page 27. We are refraining from discussing 
issues such as Solar-System tests in this paper, only concentrating on the formalism, but the 
requirements (3.26) should make it more likely that the parameter space being discussed is 
not already ruled out [103]. 

We are going to treat the expressions (3.23) and (3.24) as functions of three independ- 
ent variables: (ft, m and the expansion 9. In general one would also have to include the 
perturbations of the internal degrees of freedom of the external matter. The Einstein equa- 
tions provide some constraints between these variables. In order to obtain the perturbations 
for the hydrodynamic variables we are going to vary with respect to all these quantities. 
Eliminating the derivative w.r.t. (ft in favour of the time derivative, we obtain 

■ a 2 Q / • a 2 9\ a 2 

5S =£^ -xUO- Stf^j + m£ m A V2 , (3.27) 

5V =V^ - * (c s 2 - |) [SO - SH^j + mV m A^ , 

where A is the divergence of the acceleration defined in eq. (2.45). It expresses the spatial 
perturbation of the scalar as seen in the scalar frame (2.2), orthogonal to the perturbation 
of the proper time of that frame represented by O. We can now use (2.49) to re-express the 
59 in terms of O and A. Finally, we have to rewrite every O in terms of the momentum flux 
S, for which the definition eq. (2.21) reduces in our model (3.1) to 

H = {£ + V)Q + kA. (3.28) 

We do this since we wish to obtain evolution equations purely in terms of 5£ and we are 
able to rewrite the EMT conservation equation to retain only the combination S, eq. (2.22). 
Putting all these together, and provided £ + V ^ 0, we obtain 

„ o / 2(3 k 2 \ a 2 A faH\ 2 S k . 

6£ = e mD { 1 + £Tr^)¥ A -"{-T) h~£Tv e ' ^ 

6V ~ e \J £ + Va 2 J k 2 \ k J H £ + v[ Cs 3 

where in the coefficients we have neglected a number of terms subdominant at scales below 
the horizon. We will return to the definition of the quantity C 2 in eq. (3.34). Since we have 
no evolution equation for A, we have to eliminate the dependence of bV on it by combining 
with the equation for 5£, which gives us 

SV = CH£ - 3(c 2 - C 2 ) (f )^ l-^^-l- C>) H + (3-30) 
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where we have for the moment restored the dependence on external matter. As we have 
previously claimed in section 2.5, this is indeed a closure relation which has the structure of 
eqs. (2.48). All the parameters are purely functions of the background, apart from C 2 , which 
contains one scale as we shall see below. Nonetheless, this is a significantly more complex 
expression than that obtained for k-essence in eq. (2.46). 

We can perform the same computation for the anisotropic-stress closure relation. In the 
model (3.1), the expression for the anisotropic stress is simple 
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We can again obtain the closure relation of the correct structure by eliminating the A term 
from eq. (3.31) 
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II interpolates between a small fc-dependent quantity at large scales (effectively 0, anisotropic 
stress vanishes) and 1 at small scales, which signifies that the anisotropic stress be large. 



Let us return to the final undefined variable in (3.30), 
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where we have again neglected some terms irrelevant on subhorizon scales. This expression 
is clearly not exactly the physical sound speed of the scalar as was defined through the 
acoustic metric in eq. (3.16). However, given the assumptions (3.26), we can approximate 
the expression (3.34) as 
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(3.35) 
(3.36) 



We should stress that the above is an approximation only valid when the conditions (3.26) 
are satisfied and only on subhorizon scales. 14 We can see that as the scale k is varied, the 
quantity C 2 interpolates between two limiting values, both related to the physical sound 
speed, c 2 and c 2 — 2/3. In addition, there is a modification resulting from the Mq term 
which we should stress is defined here to be dimensionless, effectively in terms of the Hubble 
scale. 

As we will see, Mq determines one of the scales in the problem, which is equivalent to 
the Compton scale in f(R) models. It represents the degree of breaking of shift symmetry 



14 In particular, these expressions are not valid in the limit of f{R) gravity. One has to go back to the 
expression for the pressure (3.24) and recalculate assuming conditions (3.21) to be satisfied. 
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in the Lagrangian (3.1). If we make the 4> dependence of K(<f),X) and sufficiently weak, 
then only the first term of eq. (3.36) will contribute and Mq ~ C(/3/x / ). We will refer to the 
scale created by Mq in our perturbation as the Compton scale or wavelength. Interestingly, it 
appears from the definition (3.36) that M^, can be of either sign and therefore that solutions 
with increasing and decreasing Planck mass are not equivalent. One of these signs will signify 
a fast tachyonic instability for the perturbations implying that the background solution chosen 
is inappropriate. We will assume that if the effect of this mass is at all relevant then the 
background solution must be such that the mass is positive. 

In both C 2 and II there appears the same scale-dependence which changes the coeffi- 
cients in the evolution equations for the energy-density perturbations and therefore qualitat- 
ively changes their evolution. We shall discuss this scale in detail in the following sections. 

3.3.1 New Transition Scale 

Both the formula for C 2 , eq. (3.35), and for II, eq. (3.33), feature the same new transition 
scale, 

fk T \ 2 _£ + V _ 3(1 + to) 

[oHj - 2m 2 ~ 2/3 ' [6 ' 6n 
As we go across this scale, the coefficients interpolate between two values and the fluid 
transitions between two behaviours, eq. (3.33) implies that at large scales, k <C k^, the 
anisotropic stress in the fluid will be absent, while at scales k k^ we will have II ~ 1 and 
have anisotropic stress comparable to the energy density perturbation. Thus the scale k^ 
delineates the transition between what we will call the perfect and imperfect regimes. In the 
model (3.1) in the imperfect regime, the dark-energy fluid always carries anisotropic stress. 

This transition scale is determined by the relative size of the k- essence terms as com- 
pared to the second-derivative terms in the EMT for linear perturbations. The k-essence 
terms contain 5(f) and 5m, with coefficients related to the function K{(j), X). However, the 
second-derivative terms introduce terms containing (k/a) 2 5(f) and (k/a) 2 5m with a common 
coefficient k. The scale at which the two sets of terms are comparable defines (3.37). Equi- 
valently, the transition between the two regimes is defined by the scale at which the O and 
5q terms are comparable in S, eq. (2.21). Thus in the imperfect regime, it is the energy 
flow q terms that determine the evolution of the scalar-fluid perturbations. It is also worth 
stressing that it is the combination E + P that matters here, i.e. only the kinetic energy 
stored in the k-essence terms, not the potential. This means, for example, that f(R) models 
which do not have a kinetic-energy term in K(<p) are always in the imperfect regime. 

It is interesting to note that this scale is only visible through the EMT and does not 
appear in the equation of motion (3.11). The background evolution of the value of the scalar 
field <j> is determined by the background equation of motion, which has no spatial dependence. 
The evolution of the scalar 's perturbations is then determined by the perturbed equation of 
motion, which can be written as 

k 2 

5(j) + 3H(1 + -f)5(p + c 2 ^5cj) + M 2 5cp = k5p , (3.38) 

a z 

where we have neglected the gravitational potentials subdominant subhorizon and have hid- 
den the complexity of the expression in the coefficients M 2 and 7, the formulae for which we 
won't need. The gradient term for the scalar is clearly visible here and is determined by the 
physical sound speed c 2 . Provided a slow-enough variation of the coefficients, the solutions 
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to these equations are generically oscillatory. There is a change in behaviour as one crosses 
the Jeans scale, but this is the only transition. When external matter is present, it provides 
a source to eq. (3.38), which will induce an additional particular solution. This may or may 
not be dominant over the homogeneous solutions discussed, depending on the behaviour of 
all the coefficients. 

However, in order to understand how the evolution of the scalar-field perturbation 
translates into a gravitational effect which would be felt by external matter one must consider 
the relation of the scalar perturbations to the perturbed EMT (3.23). And it is here that 
a single solution of the equation of motion provides two distinct behaviours. The same 
evolution of (f> and 5(f) in time gives one behaviour to the k- essence part of the perturbed 
energy density 

5£ k 

-essence — ° "* (E^ + E m 6<j>) , (3.39) 
and a different one to the contribution from the second derivatives, 

k 2 

S£2-m - -x0-2<50 • (3.40) 

Thus the evolution of the energy-density perturbation, and therefore also the impact on the 
metric perturbations, can be completely different above and below the scale /cp, despite the 
same underlying solution for 5(f) and (f>. 

In the fluid language, we see this transition of behaviour exactly as a result of the change 
of the closure relations as we go across the scale kf. The conservation equations for the EMT 
must in the end reduce to the equation of motion for the scalar, whatever the scale we are 
looking at. Since the energy-density perturbation is a different function of 5(f) at different 
scales, the closure relations must adjust the fluid equations in such a way that the underly- 
ing solution for the scalar perturbation is always the same whatever the scale. Conversely, 
the fluid equations must predict a different evolution for the energy-density perturbation at 
different scales, given that 5(f) does not change its behaviour. This is achieved by changing 
the closure relations and the coefficients in the evolution equation for the energy-density 
perturbation (2.54), for example by introducing an anisotropic stress. 

Let us return to the discussion of the closure relations. Given the transition scale (3.37), 
we can approximate C 2 either side of it: 

C 2 ~ c 2 + M c for k <C k T , (3.41) 

C 2 ~c 2 s - 2 - + M^ ioik^>k T . (3.42) 

where we have been agnostic as to the size of the Compton term, Mq. According to our 
parameterisation, the Jeans scale is determined by C 2 + 2 /3H, eq. (2.54). Let us preview the 
implications of these approximations, leaving the detailed discussion for sections 3.3.2 and 
3.3.3: 

• If Cg S> M c , the Jeans scale lies inside the Compton scale. We can neglect M c which is 
irrelevant to the dynamics of perturbations, given an appropriate background solution. 
The Jeans term is always determined by the sound speed. 

• If c 2 <C Mq, the Compton scale lies inside the Jeans scale, eq. (3.41) implies that 
in the perfect regime the perturbation will behave as k- essence but with an effective 
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sound speed equal to Mq. The Compton mass term, as opposed to the sound speed, 
then determines the frequency of these oscillations. As we enter the imperfect regime, 
the Compton mass term will transition from playing the role of the Jeans term to that 
of a usual mass term, continuing to suppress the perturbations. There will now be a 
new scale determined by the relative size of c 2 and Mq. The physical sound speed will 
only drive the dynamics below that scale. This is effectively the picture we have in 

/(*)■ 

We have now obtained both the closure relations, eqs (3.30) and (3.32) in the necessary form 
to use in the evolution equation for energy-density perturbations (2.54) we have set up. The 
forms of parameters C 2 and II show that there are two regimes we need to consider, with 
the transition scale, &t- 

From the closure relations (3.30) and (3.32) we can now read off the values of the 
parameters defined in the general closure relations (2.48) for the case of the model (3.1). 
We have presented their values in the two regimes in Table 1. As can be seen, within each 
regime, all the parameters are indeed functions of the background quantities only, i.e. at 
most functions of time. Given this reduced set of parameters, the equivalent of eq. (2.54) for 
the model (3.1) is 15 



kr 

68 + H(8 + 3c 2 - 2wx)68 + ( C 2 + -II ) -^68+ (3.43) 

3 / a z 



+ H 2 (4(1 - IT) ^ + (3 - 2H)(5 + 3c 2 )) 58 = 0, 

Of particular note is the presence of w\ in the friction term. This term interpolates between 
and 1 and therefore reduces the friction term on scales where the anisotropic stress is present. 
As we show, this leads to slower-decaying homogeneous modes in this regime. 

3.3.2 Perfect regime, k <C k T 

In this regime, the closure relations (3.30) and (3.31) reduce to 

5V = C 2 58 + Z{C 2 -ci)^(^Pj\ (3.44) 
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C z = 4 + Mq. (3.45) 

These closure relations take the same form that generalised k- essence does (with the import- 
ant replacement of c 2 with C 2 ), see eq. (2.46). The evolution equation for energy-density 
perturbations (3.43) reduces to 

k 2 r 1 

68 + H(8 + 3cl)6S + C 2 -^68 + 3H 2 [5 + 3c 2 - 2(1 + w)\ 68 = , (3.46) 

The only correction to the standard k- essence behaviour is the fact that if the Compton 
wavelength is inside the Jeans length, M^j > c 2 , it will determine the behaviour of the 



15 We have neglected the ZU2 terms, which we can only do when \yc\ <C |1 + w\, which is implied by our 
choice (3.26). This is not the case for f(R) models (see eq. (3.22)). 
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Table 1. Closure parameters for the model (3.1) as obtained by comparing the closure relations (3.30) 
and (3.32) with the general structure proposed in eqs (2.48). The model contains two regimes, perfect 
and imperfect, separated by the transition scale fcx, eq. (3.37). In the expressions for C 2 , we have 
assumed that the external matter constitutes a negligible part of the total energy density: i.e. dark- 
energy dominates. This table is in fact exact apart from this result for C 2 , where the approximation 
(3.35) would need to be recalculated when matter is present or conditions (3.26) do not hold. In the 
presence of external matter one would need to also take into account the non-conservation of the DE 
EMT, eqs (2.52). We have provided the results for the coefficients /3j and 7.; for completeness. 



perturbations, determining the oscillation frequency of the solutions. We can rewrite the 
above as an equation for the evolution of the density contrast, 

with 5 = 5£/£. We can now obtain the standard solutions to the homogeneous equations 
in the approximation that w is slowly varying. The Jeans length in the perfect regime is 
determined by the larger of c 2 and M^, 

f k Jeans \ 2 \M 2 _ 
\ aH J C 2 

1. Above Jeans length, k <C &jeans- The tachyonic mass M 2 dominates the mass term 
allowing the DE perturbations to evolve as a power law, with modes 

6+(xa l+3w and <5_ oc a'U 1 -^ . (3.49) 

We therefore have a growing mode only for w > —1/3. For w > 1, both the modes are 
growing. In this regime, during DE domination, the fluid will only cluster when the 
universe is not accelerating. The gravitational potential is constant. 



(3.48) 



2. Below Jeans length, k 3> kj cans . At scales smaller than the Jeans length Ck/aH S> 1, 
pressure support will arrest any potential collapse, leading to oscillating solutions, 
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with Aip constants of integration and J n Bessel functions of the first kind. Note that 
the horizon is shrinking during acceleration, so the modes will eventually leave this 
regime and begin to behave as if outside the Jeans length. Depending on the hierarchy 
between c 2 and Mq, the oscillation frequency will be determined by the Compton term 
or the sound speed. 

These results match those studies in refs [64, 104, 105] for the case of k-essence. 
3.3.3 Imperfect regime, k ^> fc T 

In this regime in the absence of matter, the closure relations (3.30) and (3.31) reduce to 

£)>^-^ + <im(fy,(3,o> 




(3.51) 

It is important to stress that the closure relations (3.50) are not hydrodynamic. Firstly, the 
mass terms Mq appear in (3.50). If the Compton wavelength lies outside of the Jeans length, 
then the mass terms are always irrelevant in this regime. However, if the opposite is true, 
then the mass terms will dominate the dynamics on the outer edge of this regime. 

Secondly, it is very striking that, at small scales, where one would naively expect to 
recover the sound speed, the ratio 8V/8S can be negative when the physical sound speed c 2 
is small. If one were to take this as a definition of the speed of sound, one would interpret 
this sort of configuration as unstable. However, as we have been showing, the EMT of a 
non-minimally-coupled scalar field features anisotropic stress. This stress is very large, since 
in the imperfect regime it is of the same order as the energy-density perturbation, and it acts 
in such a way so as to precisely cancel out the term which would contribute to the would-be 
instability. 16 

This anisotropic stress from the perspective of hydrodynamics represents shear viscosity. 
This is a term appearing as a first-order correction in the gradient expansion and one would 
expect it to be small, lest it signify that the gradient expansion itself were invalid. Here, 
there are no terms higher-order in the gradient expansion to worry about: the elements of 
the EMT are are at most exactly second order in gradients. These closure relations are a full 
rearrangement of the equation of motion for the scalar field at linear order in perturbations. 

It is actually not so surprising that it is the combination of pressure and anisotropy 
perturbations that provides the support against collapse. Projecting out the spatial trace of 
the perturbed Einstein equations (the Ga components, effectively), we obtain 

/ o k 2 \ 

3> + QH$> - - (6H 2 + 4H)* = - I 5V + - -^5tt I , (3.52) 



6 The appropriate size comparison is between dimension-four quantities S£ and ^Sn, as in eq. (3.51). 
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after substituting for a scale-dependent term in Gu with the anisotropic stress. Thus the 
time evolution of the potentials reacts to the whole term in the parentheses, both pressure 
and anisotropic stress. 

We should also note that the relationship between anisotropy and the energy perturb- 
ation of eq. (3.50) is the standard result for scalar-tensor theories: the nature of the an- 
isotropic stress generated by the non-minimal coupling to gravity is such that it cancels the 
second-derivative contribution to the energy density leaving just the k-essence part of the 
energy-density perturbation. Effectively, in the imperfect regime, the dark-energy density 
perturbations in this model generate no lensing potential: 



k 2 k 2 

— ($ - tf) = 58 - -^5tt ~ 58^, 



(3.53) 



This property is an intrinsic property of this DE fluid and is always valid inside the imperfect 
regime, whatever the configuration. It is not just limited to quasi-static solutions. Some more 
general scalar-field theories have a coupling to gravity that cannot be conformally rescaled 
away and therefore the photons are sensitive to the energy perturbations of the scalar field 
(say the Lagrangians £4 and £5 in ref. [63]). 

The evolution equation (3.43) in the imperfect regime reduces to 



k 2 

58 + 3#(2 + cl)58 + c 2 ^58 + H 2 



:i + -)^ + 5 + 3c a 2 



58 = 0, 



(3.54) 



We should stress that the coefficient of the Jeans term here is just the physical sound speed 
c 2 . There are two significant differences between this equation and its equivalent in the 
perfect regime (3.46): firstly the coefficient of the friction term is reduced from H(8 + 3c 2 ) 
to H(6 + 3c 2 ) as a result of the presence of anisotropic stress. This leads to a significant 
change in the evolution of the solutions of the homogeneous modes with the scale factor 
and can affect the relative importance of the solutions to the homogeneous equation and the 
particular solutions in the presence of external matter. Secondly, in the imperfect regime, 
the effective mass of the mode has contributions from the Jeans term, the Compton term and 
the other terms depending only on the background expansion, we will therefore have three 
regimes for perturbation evolution in this case. 

We can now rewrite eq. (3.54) as an equation for the density contrast: 
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We then have three potential classes of solutions, which depend on our comparing the relative 
sizes of the terms arising from 



Compton gravitational instability 
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(3.56) 



The gravitational instability term is the imperfect-regime equivalent of the tachyonic mass 
for a perfect fluid. We have kept the standard name, even though it could be considered 
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inappropriate here since this mass is always positive, irrespective of w. As in the perfect 
regime, if c 2 3> M^j then the Compton term will have no effect on the perturbations at all. 
On the other hand, the dominance of the Compton term over the gravitational instability 
or the Jeans terms implies that we are outside the Compton wavelength of the field and 
its mass is relevant for the dynamics, shielding the interactions. At small enough scales, 
the Jeans term will end up dominating eventually, since there are always scales shorter 
than the Compton wavelength. However whether there exists a region of scales where the 
gravitational-instability terms dominate is just a question of the particular model and the 
background configuration and not scales. We should note that it is possible for the physical 
sound speed to be large enough for the Jeans length to lie already inside the perfect regime. 
In this is the case and Mq <C c 2 , the imperfect regime will only feature Jeans-term dominated 
scales. 

We can now solve for the evolution of the DE perturbations in each of the three cases: 
Jeans term, Compton term and gravitational-term domination. 

1. Gravitational instability terms dominate: 0(1 + 3w) 2 /(l + u>) 3> M^c 2 ^ 2 jk\. The DE 
perturbations evolve as a power law with 

5 + oc a( 1+3 -)/ 2 , 5_ oc a( 1+3 -) . (3.57) 

Just as in the case of the perfect regime, the homogeneous modes decay whenever the 
DE is causing acceleration, w > —1/3. However, here the mode which decays quicker 
does so only as fast as the "growing" mode in the case of the perfect fluid. There is 
now a new mode which decays even more slowly than in the case of the perfect fluid, 
c.f. eq. (3.49). This means that generically the Newtonian potential will be growing 
in this regime whenever the background is accelerating. The Poisson equation implies 
that the Newtonian potential. 

$ = 3a ^ + oca-( 1+3 "')/ 2 . (3.58) 
2 k z 

One needs to be careful, since these equations are only valid subhorizon and inside 
the imperfect regime, while the modes are leaving both and our universe is not yet 
in a state completely dominated by dark energy. However, this diametrically different 
behaviour could significantly affect late-time potentials already today. We will return 
to study these questions in detail in a separate work. 

2. Compton term dominates, M"q 3> + 3w) 2 /(l + w), c 2 s k 2 jk\ . The behaviour of the 
modes depends on the overall sign of (1 + w)Mq/ '/3. Assuming that the Compton term 
be constant, we then have two solutions, 

<ii^Mg>0, „ = U cos ^ lno + A , cos , 

a decaying oscillator or a growing mode depending on the sign. In either case the 
gravitational potential grows for an accelerating equation of state. If the Compton 
term is substantial, the tachyonic solution is a very fast instability signifying that the 
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chosen background solution itself is unstable and should not be considered a realistic 
phenomenology. 

3. Jeans term dominates, c 2 k 2 jk\ 3> Mq, (3(1 + 3w) 2 /(1 + w). The solution in this regime 
is an oscillatory function, settling to a constant amplitude when w < —1/3. The 
leading-order behaviour at late times is 

Thus contrary to the standard perfect-fluid behaviour, the density perturbation sub- 
Jeans-horizon does not decay away. Again, the Newtonian potential in this regime 
would grow while oscillating, until the mode exits the sound horizon and moves into 
regime (1) or (2). 

We have shown that at all the scales inside the imperfect regime, the perturbations evolve 
is such a way that the gravitational potential being driven by them grows. This is of course 
a frame effect: in the Einstein frame, where gravity is minimally coupled no such increase 
would be observed. However, if we assume that the observable baryons are minimally coupled 
in the Jordan frame discussed here, they would indeed see a deepening potential in this late, 
dark-energy dominated era. However, in an accelerating universe, modes are exiting the 
horizon. Thus any mode starting off inside the Jeans length in the imperfect regime, will 
exit the imperfect regime eventually and would leave the sound horizon. Thus these sort 
of solutions should be thought of as being temporary for each mode, with shorter modes 
spending longer inside the imperfect regime and therefore inducing a scale-dependence in the 
dark-energy perturbation amplitude. 

4 Discussion and conclusion 

In this paper we have proposed a unified prescription for studying the subhorizon evolution 
of linear perturbations in theories of dark energy comprising a single scalar, including, for 
example, modified gravity models such as f(R) and other models containing second- derivative 
terms in their EMTs. These EMTs are in general not of perfect-fluid form. 

k- essence is the most general scalar theory that has an EMT of perfect-fluid form. The 
evolution of linear cosmological perturbations in this class of models is already well under- 
stood. However, any further modifications to the scalar's Lagrangian, such as galileon terms 
or even a simple non-minimal coupling to gravity, generate terms including second derivatives 
of the scalar in the EMT. The presence of the second-order derivatives makes it impossible 
to express the EMT in perfect-fluid form. A standard method for dealing with such models, 
as carried out in the literature, is to turn to the equation of motion for the scalar, usually 
concentrating on the limit of negligible time derivatives for the scalar field. We have argued 
that this is in fact on one hand too specific. For example, for k-essence any two background 
solutions with the same history of equation-of-state parameter w and sound speed c 2 will give 
exactly the same evolution of linear perturbations, irrespective of the actual Lagrangian, as 
eq. (3.47) shows. One may of course prefer one of the Lagrangians over another on the 
basis of naturalness, but this is not an observable if we only have access to the cosmological 
background and the effect of linear density perturbations on gravitational potentials. On the 
other hand: the neglecting of the time derivatives could prove dangerous, if the coefficient 
of the friction term in the density-contrast evolution equation turns out to be small enough 
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to allow for modes which grow sufficiently quickly on a particular cosmological background. 
We have, for example, explicitly demonstrated that in the non- minimally coupled k- essence 
model, the DE energy-density contrast during the DE-domination era evolves in such a way 
so as to cause the gravitational potential to grow, eq. (3.57). These modes should dominate 
the gravitational potential at late-enough times on at least some of the scales. 

Our prescription allows for the study of such general models in terms of the intuitive fluid 
language instead of dealing directly with the equation of motion for the scalar field. Since 
the fluid variables such as pressure and energy density are just appropriately reorganised 
components of the EMT, they will contain all the information to which the gravitational 
field is sensitive. 

As is well known, one can obtain the evolution equation for the energy-density per- 
turbations from the perturbed conservation equations for the EMT. However, two closure 
relations need to be provided to do this: we need to relate both the pressure perturbation 
and the anisotropic-stress perturbation to the energy-density perturbation. The main result 
of this paper is the method for obtaining these closure relations, based on considering the 
actual degrees of freedom on a spatial hypersurface. 

As we have explicitly shown in our example, this method produces closure relations 
that are not those of perfect-fluid hydrodynamics. One cannot assume that the relationship 
between the pressure and energy-density perturbations is hydro dynamical, 5V = c^SS, with 
ch a hydrodynamical sound speed for pressure waves, as is frequently done. We have ex- 
plained in detail in sections 2.3 and 2.4 why even general k-essence models should not be 
thought of as a hydrodynamical fluid. The result is that the closure relations for general 
scalar-field theories take the somewhat complicated form (2.48) and cannot be reduced to a 
simple relationship. 

One may argue that this hydrodynamical relation between bV and 58 is the definition 
of the sound speed for pressure waves, ch- However, in general this is a not a physically 
meaningful quantity: in the example we have studied, where the anisotropic stress does not 
vanish, the Jeans length is not at all related to the quantity ch, but is determined by the 
speed of propagation of perturbations of the scalar field, c s , and therefore causality. As 
we explicitly show, c s is obtained by considering the effective wave operator in the scalar 
perturbed equation of motion and is not immediately apparent in the expressions for the 
pressure perturbation. Even when c s is simple and constant, the complexity of the closure 
relations we have derived shows that ch could be an arbitrary function of time and space 
and therefore not at all a useful parameter for the description of the perturbations. 

In order to aid the calculation of the closure relations, we have introduced a book-keeping 
parameterisation for them, suitable for a wide range of models containing a single scalar 
degree of freedom. This has allowed us to derive a general equation for the evolution of DE 
density perturbations in terms of these parameters which can be calculated for any particular 
class of scalar Lagrangians. For the purpose of illustration, we have shown explicitly how to 
calculate the coefficients of the closure relations (2.48) in a toy class of models, presenting the 
results in table 1. All the parameters which are non- vanishing for this class of Lagrangians 
are expressible in terms of physical quantities: w, cf, the coupling to external matter j3 
and the Compton term M^., each of which is at most a function of time but never scale. 
Non-minimally coupled k- essence was chosen as the toy model since it is the simplest one to 
contain the majority of the features which are present in more-general scalar theories such 
as galileons. This prescription can be directly applied in the study of those more-complex 
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models. 

All models containing terms with second derivatives in the EMT also contain a new 
scale &t, eq. (3.37), determined by the relative contribution of the k-essence and second- 
derivative terms in the perturbations. In our toy example, as we move across this scale, the 
DE fluid transitions from k- essence-like perfect- fluid behaviour at large subhorizon scales to 
an imperfect fluid carrying anisotropic stress at small scales. Interestingly, the scale &t does 
not appear in the perturbed equation of motion for the scalar itself, just in the EMT. On 
both sides of the transition scale, the scalar-field perturbation 6(f> evolves in the same manner, 
at least provided there isn't a scale dependence in any sources on the external-matter side. 
The changes in the pressure and anisotropy closure relations as we cross the scale k^ reflect 
simply that given the same evolution in time of 6(f), the density perturbation 6£ must evolve 
differently, since the k-essence and the second-derivative terms have a completely different 
functional dependence on 8(f). Since it is 8£ that sources the gravitational potential, this 
transition scale is physical. In addition to the transition scale, the dark energy will feature the 
known Jeans scale, controlled by the speed of propagation of small perturbations of the scalar 
field and a Compton scale related to the mass of the field. We have not studied superhorizon 
behaviour and the transition towards the horizon, but since the new phenomenology is related 
to the presence of second derivatives in the EMT, it will be irrelevant at large enough scales, 
where the standard k-essence behaviour is restored. Indeed sufficiently superhorizon the 
curvature perturbation is conversed in the standard way [59]. 

What our discussion shows is that it is more appropriate to think of the anisotropic 
stress as being related through the closure relation (2.48) to the density perturbation of the 
DE fluid rather than directly to the gravitational potentials. In general, the gravitational 
potential will also have contributions from other fluids present in the universe. The density 
perturbations in any particular fluid evolve according to the conservation equations for that 
fluid, and therefore in principle are independent dynamical variables. In our example model, 
in the imperfect regime, the closure relations imply that k 2 8Tr/a 2 ~ 6£, irrespective of the 
configuration of the external matter. It is this property of the fluid that makes the lensing 
potential independent of the perturbation in the scalar field, eq. (3.53). And it is only this 
relation that is in general independent of the other constituents of the universe. 

The only reason why in f{R) models this scale &t is not seen, is that they are a special 
subclass of the non-minimally coupled k-essence models where the k-essence terms are absent. 
Thus the DE fluid is imperfect at all scales. The Compton scale determines whether for a 
particular mode the energy-density perturbations are large or small. At scales inside the 
Compton radius they are large and therefore the anisotropic stress is large. Outside — they 
are small and so is the anisotropic stress. 

We should stress that every class of Lagrangian terms which generates second deriv- 
atives in the EMT (e.g. kinetic gravity braiding, etc.) will produce such a transition scale 
determined by the relative size of its coefficient in the background solution. Each class of 
terms will presumably cause slightly different physical effects and therefore each of these 
transitions should be potentially observable. 

Finally, we should reiterate that in order to simplify the exposition, we have mostly 
neglected in this presentation the effect of external matter. This has allowed us to concentrate 
on understanding the properties of the DE fluid itself. As a result of the non-minimal coupling 
to gravity, the DE and dark matter are coupled in this model. As we have previewed in 
the closure relations (2.48), the coupling to dark matter causes the DM perturbation to 
appear in the expression for 8V and therefore in the Jeans term. The DM perturbations 
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will in this way provide a scale-dependent source in the evolution equation for the energy- 
density perturbation 8£ , dominant in the imperfect regime. This ensures that the density 
perturbations for DE are large in the imperfect regime, providing a significant modification 
to the gravitational potentials in which dark matter propagates. 

The presence of an external energy density provides an additional complication: the way 
we defined our DE EMT means that in this case it is not conserved whenever the effective 
Planck mass evolves, i.e. in models with non-minimal coupling to gravity. As we show in the 
expressions for the conservation of the EMT (2.50) and (2.51), there are non-conservation 
terms which depend on the scalar and the dark-matter configuration. Therefore, we will 
have to provide a similar closure relation to (2.48), one for each of energy and momentum 
conservation equations. We will return to this also feature-rich discussion in a separate work. 
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A Covariant decomposition of energy-momentum tensor 

We decomposed the energy momentum tensor for dark energy according to 

= Sunlit, + V ±n V +2q^u v - ) + t^, (A.l) 
in the scalar frame defined by the vorticity-free velocity field 

V ' : " (A.2) 



We can obtain evolution equations for the fluid by taking the divergence of the EMT. It is a 
vector equation, which again can be decomposed into a time and spatial part for the observer 
moving with velocity (A.2): 

-u'V^T* = ^- + 9{£ + V) + a% + V M g* + , (A.3) 

J-A VT* = a x (£ + V) + V X V + -9q x + L x „ -f - + a x ^+ J_ A „ V ^ , (A.4) 



We have again assumed that the rotation tensor vanishes as a result of the definition of 
our velocity field. When the EMT for DE is conserved both equations (A.3) and (A.4) 
are identically zero. They then represent energy conservation (continuity equation) and 
momentum conservation (Euler equation) in the frame comoving with v?. 
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In this discussion, we are interested in analysing scalar-field theories. Since the scalar- 
field is isolated, both the energy flow and the viscous-stress tensor are formed from appro- 
priate derivatives of some scalar-field potentials, namely 

qx=±%Vpq=V x q, (A.5) 
J_»±S ~ ±^± a ^ V a V„7r = (A.6) 

V.V,, - ^ ± flu A^j 7r + + V/37T + ^$u,jy v n , 

where q denotes a potential for the energy flow while ir is the anisotropic stress, which plays 
the role of a potential for the viscous-stress tensor. They are scalar functions of the covariantly 
defined variables (f>, m (the value of "velocity" of the scalar) and potentially other scalars 
defined on the spatial hypersurface, say 6 or ^R, the intrinsic curvature of the three-slice 
(see eq. (A. 11)). The symbol A =iy- v V^Vj, is the Laplacian on the spatial hypersurface. 

Admittedly, such a simplification may not always be quite possible: the model's viscous- 
stress tensor could contain some uncontracted function of the curvature tensor C 3 )fl One 
may need to deal with these intrinsically non-scalar terms separately. The validity of such a 
decomposition should be confirmed on a model-by-model basis. 

We should note at this stage that a covariant spatial derivative acting on a function 
of only <f> vanishes (since V M /(^) = V '„(/) = —mf^ u u = 0, i.e. since the gradient 
of <f> defines the time direction). In many models, such as f(R) gravity or the non-minimally 
coupled k- essence model discussed in section 3, the anisotropic stress is only a function of (j> 
and therefore the viscous stress tensor reduces to 



' flU 



dyr(0) 
dr 



<V , (A. 7) 



and the function d7r /dT plays the role of shear viscosity. Note that depending on the direction 
of evolution of 4>, this shear viscosity can be of either sign. This is a hint that we are not 
dealing with a usual fluid. 

Assuming eqs (A.5) and (A.6), we can re-express the EMT conservation equations (A. 3) 
and (A. 4) for the case where the substitution of (A.5) and (A.6) is valid 

-vrVT^ = ^ + e(£ + V)+ 2a» V^g + Aq+ - o^V^tt , (A.8) 

±-1 vr^ = ax U + v + ^\+v x (v + ^\ + ev xq + ± x „ . (a.9) 

The gradient of viscous stress can be expressed as 

-I-Ak =~R„v ^a «" + 3*V A # + + ^V M 7r- (A.10) 



a Q V A V a vr - ^V^Avr + ^a A Avr - (3) i^V„7r , 



with the intrinsic curvature of the spatial hypersurface obtained through the Gauss-Codazzi 
equation 

(3) iV = °voP* - - l° 2 V ±- ap ±-lK • (A. u) 
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If eq. (A. 7) is satisfied then eq. (A. 10) simplifies further, 




(A.12) 



The advantage of the fully generally covariant decomposition of the EMT conserva- 
tion presented in eqs (A. 8) and (A. 9) is that it is exact under the circumstances where the 
substitution of (A. 5) and (A. 6) is valid. As we demonstrate below, obtaining the linear per- 
turbations equations for a background fiat FLRW metric is very simple. However, starting 
from these covariant equations one can also easily derive the perturbation equations for any 
other background, for example, a cosmology where the background shear is non- vanishing 
(Bianchi) or there is some sort of (potentially non-isotropic) spatial curvature. In particular, 
it is easy to see from eq. (A. 10) that for models where the anisotropic stress is more complic- 
ated than in eq. (A. 7) there is a non-trivial coupling to background spatial curvature which 
will contribute to the evolution of perturbations. 17 

We should note that the equation for the conservation of the EMT for a scalar field is 
proportional to the equation of motion for the scalar. This must be so, since the scalar is a 
single degree of freedom and providing a second evolution equation would overconstrain it. 
This means that for all scalar-field Lagrangians, eq. (A. 8) is equivalent to the equation of 
motion while eq. (A. 9) vanishes identically on all configurations. The momentum conserva- 
tion equation can in fact be shown to relate the parts of the pressure and energy density for 
the scalar, providing a sort of equation of state or, more correctly, it is the equilibrium Euler 
relation for the scalar field (see e.g. refs [45, 58]). 

B Choice of frames and boosts 

We have performed the decomposition of the EMT in what we have called the scalar frame, 
with the result that the fluid variables we've obtained (energy density, pressure, etc.) have 
the values which would have been seen by an observer moving with velocity defined in eq. 
(2.2). A different observer will in general observe a different energy density, etc. (see also 
the discussion in ref. [107]). 

The k- essence class of models has EMTs of perfect-fluid form. The energy- flow q^ in 
the scalar frame vanishes and thus the scalar frame is also the Landau-Lifshitz frame (or rest 
frame). This makes the choice of frame for the analysis very simple: the scalar frame is the 
natural one. However, in general scalar- field models (for example, see ref. [58] or section 3), 
the energy flow in the scalar frame is non- vanishing. In these circumstances, the LL frame 
is distinct from the scalar frame and one needs to choose. Let us compare the physics in the 
two frames. 

The two frames are related by a (space-dependent) Lorentz boost. However, in general 
there isn't a simple transformation between the two frames. For the purposes of cosmology, 
q^ and 77j„ in the background vanish because of the symmetry of the FLRW metric, so we 
are able to make a perturbative statement about this transformation. 

The Lorentz boost we require is a transformation from the orthogonal vector pair (u^, g M ) 
to the orthogonal pair (U^, Q^) defined by the linear transformation 




(B.l) 



17 



For example, if the action for the scalar contains a term G^V^V^ [106]. 



-43- 



where cp 1 is the unit vector in the direction and the rapidity a is defined as 

« 2 = f% , (B.2) 
£ + 7? v ; 



~cT^F, , (B.3) 



or 

£ + P 

to first order in a. The frame comoving with velocity {7 M is now the LL ("rest") frame. In 
particular, the EMT decomposed in this new frame to leading order in a can be written 




with i-ui/— + U^JJ V . This implies that up to first order in the rapidity a, an observer 
comoving with velocity will observe the same energy density, pressure and viscous stress 
as one comoving with u M : we have 

Just as required, the energy flow in this frame vanishes. However, the decomposition (B.4) 
shows that the observed LL fluid variables are in fact modified with respect to the scalar 
frame already at the second order in a. 

One should note at this point that the above is not a gauge-dependent statement. 
Frequently the question of gauge and frame choice are conflated in cosmology. They are 
in principle very different. Gaug GS £11*6 cl choice of the form of the metric perturbations 
which are all equivalent as a result of diffeomorphism invariance of general relativity. They 
do not impact observables. Since the Newtonian potentials $ and are gauge invariant 
and are equal to the potentials <pL and ipL in the longitudinal gauge, one should consider the 
perturbed equations (2.18) and (2.19) to be gauge-invariant themselves, provided that each of 
the quantities is interpreted as the gauge-invariant variable with appropriate transformations 
given different gauge choices. The choice of frame, on the other hand, is a statement about 
observables. As can clearly be seen in eq. (B.4), the quantity that the observer moving with 
will call the energy density is a scalar, but it is a different scalar to the energy density 
seen by u^. 

Turning to first-order cosmological perturbation theory, we can find the spatial velocity 
divergence for corresponding to (2.20), 

©ll = ikiSU* - ®x ~ o k l 5q ^ + iJ-^Qx , (B.6) 

a 2 (£ + T) £ + V K ' 

and therefore we can rewrite the conservation equations (2.18) and (2.19) in the Landau- 
Lifshitz frame, 

5£ LL + 3H (5£ LL + 5V LL ) + {£ + V) (e LL + 3$) = , (B.7) 

cr 3 a q 



{£ + V)[ 9 LL + 2if9 LL - -J* - -^5V LL + QllV - ^5tt ll = . 
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where the subscript LL on the fluid variables signifies that these are the quantities observed in 
the LL frame and just happen to be to equal to the scalar-frame equivalents up to first order 
in a. Unsurprisingly all that has happened is that the terms containing q have disappeared: 
the form of these equations must be invariant under Lorentz boosts. 

One could at this stage walk away having decided that one should always perform 
cosmological perturbation theory in the Landau-Lifshitz (rest) frame since it appears com- 
putationally simpler as a result of having eliminated the additional variables q. However, 
this is dangerous for two reasons: 

• The velocity field is now no longer vorticity-free. The scalar-frame velocity field 
was a (normalised) derivative of a scalar field and thus the antisymmetric rotation 
(twist) tensor vanished. This is no longer true for C/ M and therefore the scalar-field 
fluid can carry vector perturbations in this frame — albeit only at higher orders in 
perturbations. One may not be particularly worried about this since at higher orders 
all the perturbation types mix anyway, but it does provide additional complication, 
especially with regard to choosing a good Cauchy surface. One would also need to add 
the rotation tensor terms to eqs (A. 3) and (A. 4) and would not be able to perform the 
decomposition leading to (A. 8) and (A. 9). 

• Possibly more importantly, eq. (B.2) shows that the rapidity a required to boost to 
the LL frame diverges when the equation of state of the fluid approaches w = —1. The 
physical origin of this is that the vanishing enthalpy of the vacuum cannot compensate 
for the energy flow q^. This means that as the fluid approaches the vacuum equation 
of state, the perturbative relations between the fluid variables in the LL and any other 
frame breaks down: in the LL frame, the energy-density functional form itself contains 
corrections (the leading one can be seen in (B.4)), which are order one close to w = — 1. 
This also means that order-one vorticity would be observed in the LL frame. 

The last point implies that if there is any danger of having an equation of state close to that 
of vacuum, one should not decompose the EMT in the LL frame and attempt to evolve these 
hydrodynamic variables. On the other hand, any other choice of frame is in fact equivalent 
up to first order in the difference in velocities and for the purpose of linear perturbation 
theory one can directly use results from one frame and apply them to others. 

One could of course claim that such a setup is not possible: after all one calculates the 
hydro dynamical EMT by averaging quantities in the rest frame of the particles making up 
the fluid in the first place. Why should it be possible to produce an EMT which one could 
not put back in the rest-frame? As we argue in section 2.4, general scalar-field theories do 
not have EMTs which behave as in hydrodynamics and indeed have been shown to be able 
to produce a fluid vacuum configuration than nonetheless carries momentum in the spatial 
direction (e.g. ref. [58]) 

It is worth noting here that the practice assigning the meaning of e.g. the energy density 
to components of the EMT can also be understood in the framework of the decomposition 
(2.14). First, one should realise that it is impossible to find a frame where the definition of 
energy density is exactly —T° given an arbitrary metric. The best we can do is to pick the 
time direction to be that of the coordinate time of the background, i.e. 
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where the square root of the metric is necessary for the appropriate normalisation of the 
velocity. In this frame we can then obtain the following decomposition of the EMT: 

-g™T 00 = -T + g 0l T i0 , (B.9) 
5{^°°(T oi + g 0i £) , 

$1% (Tij - V gij - g 00 (£ - V) g oi g oj - g 00 (T 0l g 0j + T 0j g 0t )) . 

This then matches the components approach, but only provided that the (Oi) metric elements 
vanish. Even in this case there is an additional metric factor multiplying q\. However, 
given scalar perturbation of FLRW in the Newtonian gauge and neglecting second-order 
corrections, this decomposition yields the same energy density, pressure and viscous stress as 
the decomposition in the scalar frame. 

Another aspect worth mentioning is what happens in the case of an external EMT. 
For example, if we consider the scalar field to model dark energy, dark matter would be an 
additional external perfect fluid. The external EMT will contain its own rest frame with 
velocity, say, 

= Pw± v H v u + Pext i-^ , -Lj^= V + v ^ v v ■ (B.10) 

However, we have to pick just one frame in which to do the decomposition, and it is convenient 
to make that the scalar frame. This means that the fluid variables describing the dark matter 
are not those observed in its own rest frame, but appropriately transformed quantities. For 
example, the external pressure observed in the frame would then be 

i&t = ^ext (2 + (u ■ v) 2 ) - ^ cxt (l - (u ■ v) 2 ) . (B.ll) 

The various properties of the external perfect fluid are defined in its own rest frame (just 
as they are for k- essence, for example). Thus in general one needs to keep track of how 
the external matter appears in the frame in which we are decomposing and perform the 
perturbations. For example, eq. (B.ll) show that external dust, which has an equation of 
state w ex t = in its own frame, would have a non- vanishing pressure in the scalar frame. 

However, just as in the case of the perturbations of the scalar, if both the fluids share 
the FLRW background, the corrections to energy density, pressure and viscous stress arising 
from the difference in the rest-frames appear only at second order, thus for example the 
energy density perturbations 5p ex t in the external rest-frame is going to be the same as in 
the scalar frame and we can drop the superscript U used in eq. (B.ll) to differentiate the 
two. At linear order the only difference will be in the energy flow: since the two fluids do 
not share the same rest frame, the external matter will have a non-vanishing velocity from 
the point of view of the scalar: 

Sq^t = T ex V ± ^ + pcxt) ^ {8v . _ 5ui) _ (B 12) 

or if we take the divergence, 

iki$ql xt = (Pext + Pexp) (©ext - ©) , (B.13) 

where O ex t = di5v l is the analogously defined spatial divergence of the external matter 
velocity perturbation. 
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To conclude: the choice of frame is in general important when one is worried about 
higher orders in cosmological perturbation theory. If we only concern ourselves with linear 
perturbations, then any frame where the observers while unperturbed are comoving with the 
FLRW background observers will give the same results. However, for scalar-field models with 
second-derivatives in the EMT, the Landau-Lifshitz (rest) frame is singular when w = —1 
and perturbation theory loses validity as the fluid approaches the vacuum-energy equation 
of state. Thus, somewhat counter intuitively to expectations, one should use any frame but 
the rest frame to study linear perturbations in cosmology. 
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